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Abstract 



The BEC-BCS crossover problem has been intensively studied both theoretically and experimentally largely 
thanks to Feshbach resonances which allow us to tune the effective interaction between alkali atoms. In 
a Feshbach resonance, the effective s-wave scattering length grows when one moves toward the resonance 
point, and eventually diverges at this point. There is one characteristic energy scale, 6c, defined as, in the 
negative side of the resonance point, the detuning energy at which the weight of the bound state shifts from 
prcdominatedly in the open-channel to predominated in the closed-channel. When the many-body energy 
scale (e.g. the Fermi energy, Ep) is larger than Sc, the closed-channel weight is significant and has to be 
included in the many-body theory. Furthermore, when two channels share a hyperfine species, the Pauli 
exclusion between fermions from two channels also needs to be taken into consideration in the many-body 
theory. 

The current thesis addresses the above problem in detail. A set of gap equations and number equations 
are derived at the mean-field level. The fermionic and bosonic excitation spectra are then derived. Assuming 
that the uncoupled bound-state of the closed-channel in resonance is much smaller than the inter-particle 
distance, as well as the s-wave scattering length, a^, we find that the basic equations in the single-channel 

crossover model are still valid. The correction first comes from the existing of the finite chemical potential and 
additional counting complication due to the closed-channel. These two corrections need to be included into 
the mean-field equations, i.e. the gap equations and the number equations, and be solved self-consistently. 
Then the correction due to the inter-channel Pauli exclusion is in the order of the ratio of the Fermi energy 
and the Zeeman energy difference between two channels, Ep/rj, which can be analyzed perturbatively over 
the previous corrections. 

Fermionic and bosonic excitation modes are studied. Similarly as the mean- field result, the basic structure 
follows that of the single-channel model, and the correction due to the inter-channel Pauli exclusion can be 
treated perturbatively with expansion parameter in the order of Ep/r]. In the bosonic excitation, a new 
out-of-sync phase mode emerges for the two-component order parameters. It is nevertheless gapped at the 
the pair-breaking energy. 
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Chapter 1 

Introduction 



From a methodological view, a physical system would be very desirable for developing and verifying a theory 
if it can be described with as few parameters as possible, and each parameter is as tunable as possible. One 
such system is the ultracold dilute fermionic alkali gas with the Feshbach resonance. Dilute fermionic alkali 
gas was cooled into degenerate region in 1999 [DeMarco and Jin, 1999| . Not long after that, superfluidity 
was observed for such systems in 2003 [Regal et al., 2003 . In dilute ultracold fermionic alkali gas, it is 
sufficient in many phenomena to describe an atom-atom interaction with one single parameter, namely the 
s-wave scattering length, as, because the gas is very dilute and experiments are performed at very low 
temperatures. 

The other desirable property is the ability to tune the effective interaction strength, i.e., in the present 
case, the s-wave scattering length, through Feshbach resonances. One atomic energy level splits into 
several hyperfine levels under a magnetic field, due to the hyperfine interaction between nuclear spins and 
electronic spins. Hyperfine spin indices provide a good set of quantum numbers for a single atom. In the 
theory of atom-atom interactions, a channel refers to a specific configuration of hyperfine spin indices of one 
atom pair. Different channels usually have different magnetic moments and therefore have different Zeeman 
energies in the presence of a magnetic field. The difference of the Zeeman energies can then be tuned by 
changing the magnetic field. In addition, a channel is no longer an exact eigenstate when the atom-atom 
interaction is taken into consideration because the interaction is mostly due to the overlap between electronic 
parts of the two-atoms' wave- functions. As a result, channels are hybridized. The effective potential of each 
channel is also different. The potential of one channel may be deep enough to sustain a bound state. 
This channel is then called "closed-channel" . For a certain magnetic field, the energy level of this bound 
state might be close to the zero-energy threshold of the other channel, usually called open- channel, and 
two channels thus strongly hybridized. The low-energy scattering properties in the open-channel is then 
dramatically modified. In such a situation, two atoms approaching each other in the open-channel may 
"spend a certain amount of time in the closed-channel" and then reemerge in the open-channel. Atoms in 
the open-channel seem to feel an enhanced effective interaction. This phenomenon is known as Feshbach 
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resonances. We will present a quantitative analysis about alkali gas in Chpater [5] and Feshbach resonances 
in two-body context in Chapter [31 

A very desirable property of the Feshbach resonance is that the effective interaction is tunable experimen- 
tally through the Zeeman energy difference between channels which is in turn tunable through instruments 
such as a magnetic field. This unique ability gives physicists a rare opportunity to study a many-body sys- 
tem under various interaction strengths, and thus connect different physics originally developed separately. 
Particularly for the fermionic gas, there are a series of theoretical works about uniform treatment of BEC 
and BCS since the 1960s [Eagles, 1969l|Leggett, 1980l|Nozieres and Schmitt-Rink, 1985[|Randeria,"T995] , for 
which dilute ultracold fermionic alkali gas with Feshbach resonances provides the perfect testing grounds. 
Indeed, these theories works quite well qualitatively. 

> 0, 5 < / 

a, < 0, (5 > 



zero energy scattering state / ttg ^ OO 





closed-channel dominoed open-channel dominated 




bound state 



Figure 1.1: Energy levels in a Feshbach resonance 
S is the energy detuning from the resonance point, where the resonant point is 
defined as the point where the open-channel efi'ective s-wave scattering length di- 
verges, fls — ^ ±oo. The horizontal line stands for the zero energy s-wave scattering 

state, ip ^ —, which exists for any detuning. The lower line stands for the 

real bound state, which only exists for negative detuning {S < 0, > 0). The 
dash hne stands for the (uncoupled) closed-channel bound state. An interesting 
point to notice is that the real bound state appears earlier than the cross point 
of the (uncoupled) closed-channel bound state level and zero energy. Another im- 
portant point to notice is the negative detuning —Sc- When the negative detuning 
is smaller than Sc, this real bound state is composed mostly with atoms in the 
open-channel and vice verse. See Chapter [3] for details about K, and Sc. 



The two-body theory of the Feshbach resonance has a characteristic parameter, Sc, defined as the detuning 
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energy at which the weight of the bound state shifts from predominated in the open-channel to predominated 
in the closed-channel (see Fig. II. ip . Naively speaking, in the negative detuning side of any resonance (i.e. 
(5 < 0), particles should mostly stay in a (virtual) bound-state of the closed-channel (or "virtual state" in 
some other resonances). However, at the resonance point of the Feshbach resonance (os — ^ ±00), atoms 
are mostly still in the open-channel, and they do so down to a negative detuning S ^ ~Sc- Only when the 
detuning from resonance is much far away than Sc, do atoms have the majority weight in the closed-channel. 

When considering a many-body problem, an important question is how this energy scale, Sc, compares to 
a typical many-body energy scale, namely, the Fermi energy of the free fermionic atoms. In the region not 
too far away from resonance {\S\ <^ Sc), the closed-channel weight is negligible if the Fermi energy is much 
smaller than 6c, (i.e., broad resonance). Crossover experiments are usually performed at detuning not too far 
from the resonance, and hence the closed-channel can be safely ignored at the many-body level. Eventually, 
when the detuning is too far away, \S\ ^ Sc, the bound state is almost like the uncoupled closed-channel 
bound state with a little dressing from the open-channel. We nevertheless do not concern such cases for the 
broad resonance because crossover phenomena have already be well covered in both BEC and BCS ends with 
\S\ <^ Sc- The problem can be well-described as a two-species fermion system with a tunable interaction. 
The Feshbach resonance indeed serves as a simple "magic" knob to change the interaction strength. The 
original theories developed on single-channel models apply to this case directly. This is also the situation for 
two popular experiment cases (^Li atoms at 834G, ''"K atoms at 224G). Many theoretical works have been 
developed using the single-channel model as these original works or using the tow-channel model within the 
broad resonance assumption (e.g. Holland et al., 2001 |Ho, 2004|Fuchs et al., 2004| ). On the contrary, when 
the Fermi energy of the free fermionic atoms is comparable to or even larger than Sc, the closed-channel has 
to be included at the many-body level even for small detuning. Such a situation, previously considered in 



some works Gurarie and Radzihovsky, 2007 , is the focus of the current thesis. 

Nevertheless, one crucial simplification comes from the fact that relevant uncoupled closed-channel bound 
state is tightly bound, with spatial extension much smaller than many-body scales, e.g. the interparticle 
distance, (but often larger than the potential range). This fact enable us to treat the Pauli exclusion 
between two channels perturbatively. It is not necessary to handle all the fermion species simultaneously, 
which probably requires quite different techniques other than those discussed in this thesis. 

To complicate the problem even further, real experimental configurations often have one common hy- 
perfine species between the two channels. There are three hyperfine species in the two channels instead of 
four species (two for each channel). Two most common setups (^Li at 834G, ''"K at 224G) both contain 
three species of fermions although they are the broad resonance. The Pauli exclusion principle prevents 
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atoms of both channels from occupying the same level simultaneously because of this common species. The 
inter-channel Pauli exclusion has no counterpart in two-body physics. This peculiar effect in many-body 
crossover problems has received little theoretical attention up to now. Nevertheless, narrow resonances do 
exist [Chin et al., 2010] and it is not inconceivable to perform many-body experiments using such resonances. 
The central concern of this thesis is about these situations. 

Roughly speaking, turning from two-body systems to many-body systems brings three effects into the 
original two-body problem. The first effect is closely associated with the Fermi energy: For a many-body 
fermionic system at low temperature, most fermions are inactive; only the fermions close to the Fermi 
surface participate in the interaction processes. Therefore, energy often needs to be measured from the 
Fermi surface instead of from zero as in a two-body situation. This aspect has been extensively studied 
previously [Gurarie and Radzihovsky, 2007| . 

The second effect is about counting. Unlike in the single-channel problem, there are two relevant densities 
in the two-channel problem: the density of atoms in the open-channel, tIo, and the density of atoms in the 
closed-channel, ric- When the closed-channel weight is small (broad resonance), it is all right to treat 
the total density as the same as the open-channel density. However, in the narrow resonance, where the 
closed-channel weight is not negligible, counting becomes complicated. Extra care is required to specify 
which channel the quantities, such as "density", belong to. This aspect has been also extensively studied 
previously [Gurarie and Radzihovsky, 2007| . 

The last effect is unique for the three-species problem, where one common species is shared by both 
channels. The phase spaces of two channels are overlapped because of the Pauli exclusion caused by the 
common species. This effect is controlled by the wave-function overlapping of states in the two channels. 
A rough estimate of this overlapping can be made: The uncoupled closed-channel bound-state which is 
in resonance with the open-channel zero energy threshold has relatively small spatial extension, Oc. Its 
binding energy Ef, is close to the Zeeman energy difference between two channels, 77. On the other hand, 
fermions in the open-channel fill the lowest momentum states up to typically the Fermi energy, Ep. By a 
simple dimensional argument, the ratio Ep/rj must control the overlapping effect. This effect has not been 
addressed in any theoretical work to this author's knowledge. How it modifies the many-body picture is the 
central topic of this thesis. 

We can see this from a slightly alternative aspect using two-fermion molecule gas for the uncoupled 
closed-channel bound state. We assume that the molecule size is Qc and the total number of molecules is 
N. Assuming further that the bound-state is close to threshold, the bound-state wave function can then be 
written as A/{k'^ + k^), where h?K^/2m — Eb, (see Appendix IB. 3p . The prefactor "A" can be determined 
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by normalization, 'Yl!k=o 1^1^ ^ ^ ■ Now we consider all atoms in a typical many-body scale, e.g. the Fermi 
energy, Ep^ which is going to overlap with levels occupied in the open-channel. Usually, the Fermi energy 
is much smaller than the energy scale of the closed-channel bound state, Ep <C E\j. The total number of 
atoms in [0, Ep] is roughly N ■ [kpacY , which is much smaller than N . This means that in the two-channel 
problem, the low momenta, (fc < kp)^ are still dominated by the open-channel component even when the 
total number of atoms in the closed-channel is comparable or higher than the total number of atoms in the 
open-channel because atoms in the closed-channel are mostly in high-momentum states. 

The present thesis is divided as follows: Chapter [2] to S] review several important concepts used in the 
thesis. Chapters[S]and[n]then present my main work and Appendix \K\ lists an earlier attempt using a roughly 
equivalent but less-flexible approach. 

More specifically, Chapter [5] briefly reviews dilute ultracold alkali gas. Section [^31 in particular examines 
the idea of "universality" , which is one of the central ideas in our treatment of the two-channel model. 
Chapter [3] goes over the Feshbach resonance in two-body physics and the concept of the narrow (broad) 
resonance is introduced. Chapter 2] reviews the single-channel BEC-BCS crossover problem as well as the 
path-integral approach solving it. This chapter serves as the starting point for the solution of the two- 
channel model. After these reviews. Chapter [S] and Chapter |5] present my work on the three-species narrow 
Feshbach resonance within a many-body path-integral framework, in detail. Chapter [5] discusses the mean 
field result while Chapter [S] discusses fermionic and bosonic excitations. An earlier attempt based on the 
BCS ansatz approach in mean-field level is given in Appendix [Xj Chapter [7] discusses our procedures and 
their conclusions. 
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Chapter 2 

Dilute ultracold alkali gas 



Since the 1990s, dilute alkali gas has been cooled into quantum degenerate region where the thermal de 
Broglie wavelength (--^====) is comparable or larger than the interparticle distance. Not long after suc- 
cessfully cooling the bosonic atoms, fermionic alkali gas was also available in the degenerate region. Because 
of the ultra-low temperature (in the order of nK), and the extreme diluteness (10^^ ^ 10^'^cm^'^), the atoms 
are mostly free except when they are close. This particular property simplifies theoretical analyses tremen- 
dously (see Sec. 12.31 for details). In this chapter, we review a few aspects of the dilute ultracold alkali gas 
that are closely related to the current thesis. 

2.1 A single atom and its hyperfine levels 

In experiments on ultracold alkali gas, a magnetic field (B) is the most common physical quantity to 
manipulate. Let us first study an isolated atom in the presence of a magnetic field. An alkali atom has 
only one electron in its outer shell. All the rest electrons are in the filled inner shells which has no total 
magnetic moment. So we only need to consider the spin of the outermost electron, S, for interaction with 
the magnetic field. In addition, a magnetic field also interacts with the atom's nuclear spin, I. The full 
spin-part Hamiltonian is 

Hspin = AI • S — /igB • S — /i„B • I 

(2.1) 

= AI • S - HeBSz - fJ-nBIz 

The first term with a characteristic energy A describes the hyperfine interaction, while the next two terms 
describe Zeeman energies of the outer electron and nuclei respectively, /ig is the electronic magnetic moment, 
while /i„ is the nuclear magnetic moment. In the second line, we take the direction of the magnetic field as 
the z-direction. This Hamiltonian can be diagonalized by introducing the total spin 

F = S-HI (2.2) 
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When the magnetic field is zero, the two Zeeman energy terms in the above Hamiltonian vanish. Tlius, {F, Fz) 
are good quantum numbers and all states with the same F are degenerate. When the magnetic field is finite, 
{F, Fz) cannot diagonalize the Zeeman energy terms, and therefore are no longer good quantum numbers. 
Nevertheless, we can still label the atomic states with these two numbers via the adiabatic connection to the 
levels in the zero magnetic field. For a finite magnetic field, except for states with the highest and lowest 
Fz, namely, ±F, each other state is a mixture of different {Sz,Iz) or {F,Fz). Fortunately, S is just equal to 
1/2 for an alkali atom. So each state is mixed with at most two sets of quantum numbers {Sz, Iz)- At high 
magnetic fields, the first hyperfine coupling term in Eq. (j2.ip is dominated by the last two terms of Zeeman 
energies and the eigenstates are approximately described by the quantum numbers (S'2,/2). (See Fig. 12.11 ) 



E n);hlMHzL 




Figure 2.1: Hyperfine structure of a single ^Li atom 
Levels are marked with F and Fz (see Footnote [T] in page [7]) 



2.2 Two-body interactions 

Things would be very boring if there was only the single-atom Hamiltonian. Before discussing the interac- 
tion itself, let us introduce one important concept related to it. The term "channel" is used to refer one 
configuration of hyperfine spins for one atom pair in interaction, 
are good basis for non-interaction pairs. A pair of atoms in one channel would stay in this channel forever 
in the absence of interactions between atoms. Considering the present case of alkali atoms, two alkali atoms 



L lie lerm cnannei is useu uo reier one 
F'^'\F^^'^) . Q Channels 



^Remember that {F,Fz) are only labels; they do not stand for the total angular momentum unless there is no magnetic 
field. 
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mostly interact through the overlapping of their electron wave-functions in the dilute limit. Thus, besides 
the atom-atom distance, the interaction is mostly a function of electronic spins, with negligible dependence 
on nuclear spins. Schematically the interaction can be written as 

V^f{r)+g{r)Si-S2 (2.3) 

Here the spatial part of the interaction is coupled with electronic spins. The hyperfine levels that diagonalize 
the single atom Hamiltonian are no longer eigenstates for this interaction. In another word, the interaction 
has non-diagonal terms between channels and therefore hybridizes them. Instead, states with definite elec- 
tronic spins form a good approximate basis for the atom-atom interaction. Nonetheless, most experiments 
are performed in the so-called high-field region, where the Zeeman energy terms dominate the hyperfine 
interaction term in the single-atom Hamiltonian (Eq. 12. Recalling that one hyperfine state is at most 
mixture of the two {Sz,Iz) states. In the high-field, one of them always dominate the other; therefore the 
electronic spin Sz is approximately a good quantum number for a hyperfine level, and the original hyperfine 
levels (channels) serve as a good starting point as the zeroth order approximation. When the hybridization is 
taken into account, multichannel scatterings are possible. The most interesting thing in these multi-channel 
scatterings is the possibility of a resonance. One of them is the so-called Feshbach resonance: The potential 
in one channel may be deep enough to sustain a bound-state. This channel is named "closed-channel" . When 
this bound-state energy is close to the zero-energy threshold of another channel, named "open-channel" , the 
low-energy scattering properties in the open-channel are dramatically modified. This resonance turns out 
to be extremely useful in cold atom experiments. Chapter |3] reviews its theory for a two-body system. The 
more involving many-body problem is then the central theme of this thesis and is discussed in Chapter [S] 
and [5] in detail. 

Let us illustrate above discussion in one example. In a common experimental setup for ^Li (Fig. 12. ip . 
experiments are usually prepared with atoms in the two lowest hyperfine levels: described by the direct 
product of two states, \F — ^,Fz — —\) ®\F ~ ^,Fz ^ +\)- This is a good approximation until the two 
atoms are very close. Recalling that the atom-atom interaction (Eq. 12. 3p conserves the z-component of the 
total angular momentum, F^'^'^ -f Fi'^\ this channel mixes with four other possible channels of the same total 
z-direction angular momentum, i.e. F^^^ + F^'^^ = 0: | |i ® 1 5' 1 1' +1) ® 1 1' 

1 1, +^)® I f , — ^) (All states are labeled as |F, F^)). Various resonances can take place. Note that close to the 
resonance, it is normally sufficient to consider only the channel that is in resonance, while neglect all other 
channels. Another important aspect is whether the two channels share one single hyperfine species or not. In 
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the former case, totally three hyperfine species are involved while in the later case, four hyperfine species are 
involved. The closed channel in the most studied resonance with a magnetic filed close 834G, is approximately 
If,— 5) <8) and the resonance is a three-species resonance [Zhang, 2009[|Chin et al., 2010| . 

2.3 Universality, Bethe-Peierls boundary conditions, the s-wave 
scattering length, and two-body density matrices 

One important aspect of the interaction in dilute ultracold alkali gas is that for many purposes, it is sufficient 
to characterize the interaction by a single two-body parameter, namely, the s-wave scattering length, Os, 
because both the density and the temperature, are very low. This is often interpreted as we can replace the 
real potential with a pseudo potential, U{r) = ^^^j^5{r) [Pethick and Smith, 200T||Leggett, 2001| . Never- 



theless, an alternative interpretation of [Leggett, 2001[[Tan, 2008a| Tan, 2008b Combescot et al., 2009] is 



more useful in this work. For a short-range potential, where the potential range, Tc, is much smaller than the 
average interparticle distance, ao, particles are free-like in the majority of the time. They only interact when 
two particles are close to each other. We can thus schematically divide the space into two domains: V, where 
any two particles are more than away from each other; and otherwise, I. Most physical quantities would 
be very easy to calculate if only considering the free part, V. In a low-energy (ultracold) dilute system, we 
only need to consider pair-wise interaction while neglect all the three-body or more-body interactions. In 
addition, V takes the majority of the space due to the same reason. The effect of the potential on wave- 
function in the short-range region, I, can be taken simply as a boundary condition on the wave-function 
in the free part 2?, ^p{r) i'oi'f)- For an isometric ^po{r), the lowest order in the radial coordinate r is 



^. y Including the next order, a constant, we have ipoi''") 7 ~ ^ barring the normalization. All these 
consideration gives us the simplest non-trivial boundary condition on a wave function 

^(r)^A(^-^) (2.4) 



which is also known as the Bethe-Peierls boundary condition [Bethe and Peierls, 1935 . is fully determined 



by two-body physics. This simple boundary condition applies to two-body, few-body, as well as many-body 
systems, and has been proved to be a very powerful tool in solving various problems. 

Eq. 12.41 coincides the zero-energy s-wave scattering wave function, which explains the name of parameter 
"os", the s-wave scattering length. Nevertheless, we did not mention anything about zero energy so far. 



^The extra - factor is there for radial wave function in 3D. 



9 



although Qs is defined for zero energy in the scattering theory context. In fact, this boundary condition (Eq. 
I2.4[) apphes generahy to any low (positive or negative) energy solutions as long as the energy involved is much 
lower than the energy scale in the interaction domain I. Hence, this boundary condition applys to close-to- 
threshold bound states as well. The s-wave wave function of a weak bound-state reads ij{r) — ^e^''/''^ in 2?, 
which matches the Bethe-Peierls boundary condition with a positive as for r ^ a^. The exponential decay 
factor of the wave function, a^, is directly related to the binding energy. Eh, with the often cited relation. 

- (2.5) 

Here is the reduced mass for center of mass, which is equal to half of the atom mass for a pair of the 
same atoms. This immediately clears one often confusing and counter-intuitive fact, that a positive is 
associated to a bound state. In the standard scattering theory, a positive as is usually associated with a 
repulsive interaction, which obviously does not support a bound stated 

Eq. 12. 4[ does not fix the normalization on the wave function. This normalization factor, encapsulating 
effects from particles outside the immediate interacting pair, appears in many physical quantities. In a 
dilute and low-energy system, its square is proportional to the so-called integrated contact intensity, C , 
introduced by Tan |Tan, 2008a"l|Tan, 2008b||Combescot et al, 2009] . For the hmit 0, the integrated 



contact intensity, C, and the s-wave scattering length, as, are sufficient to describe several important physical 
quantities, such as internal energy. A particular useful one for this thesis is the limit at high-momentum of 
the particle number distribution of particles, 

Note that here high-momentum does not mean the absolutely high-momentum, it means lower than the 
characteristic momentum of potential l/rc but higher than any other scale, 1/ao,... Indeed, when the 
short-range approximation and low-energy approximation apply, we expect that the two-body correlation 
at high-momentum (3> kp) does not change much from a two-body system to a many-body system. In this 
high-momentum region, we can always use the two-body wave function as good approximation. 

In many-body physics, various physical observable quantities are related to one set of quantities, namely, 
density matrices, (^'^(xi) • • • ^'^(xAr)^'(a;'^) • • • \I'(a;']^)). In fermionic system, the one-body density matrix 



^This seeming paradox can be resolved carefully within scattering theory as follows. In the scattering theory, only when the 
interaction is weak, and phase shift as well as as are small, we have the fact that a repulsive interaction leads to a positive phase 
shift and a positive ; while an attractive interaction leads to a negative phase shift and a negative ■ No simple relationship 
of signs holds for a strong interaction, where a bound state might form. In fact, the phase shift changes as much as 2it when a 
bound state starts to form; therefore, is large and changes sign over the threshold. 
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(A^ = 1) for an interacting system is often very close to the one for the free particles, although the difference 
can be important for some theories, such as Landau Fermi liquid theory. A two-body density matrix {N — 2) 
is often used because it is often different from the two-body density matrices of the free fermions qualitatively, 
such as in the case of BCS pairing. Formally, we can decompose a two-body density matrix into an orthogonal 
basis 

(Vl/t(xi)vl/t(:^2)vl,(y2)vl/(yi)) =^C„4(xi,X2)0„(yi,y2) (2.7) 

n 

When one or a few C„ are macroscopic, some special quantum phenomena often emerges. Especially when 
only one parameter, Co, is macroscopic, the system can often be interpreted as one macroscopic wave 
function, (j>o{xi,X2), (which is often called order parameters). [Leggett, 2006| This can serve as a starting 
point for several phenomena, such as BCS superconducting,... 

Zhang and Leggett developed independently another universality theory based on two-body density 
matrices Zhang and Leggett, 2009] , which actually take a more general form of boundary conditions than 
the Bethe-Peierls boundary condition with as- They asserted that for a short-range potential and a low 
temperature, as in the case of dilute ultracold alkali gas, the basis wave functions 0„ in Eq. (|2.7p follows 
the two-body wave function at short-range. This is actually similar to Eq. (|2.4p . Instead of requiring the 
simplest form of ipo given in Eq. (|2.4[) . they require a more general wave- function that solves the Hamiltonian 
at two-body level. Not surprisingly, many physical properties are determined by the normalization factors 
as using the Bethe-Peierls boundary condition. 

In this thesis, similar idea is used. We assume that the closed-channel correlations follow the its two- 
body bound-state wave- function at high energy (i.e. short-range). An open-channel correlation follows its 
two-body wave-function at short-range, but its intermediate range does not do so because of the sensitive 
nature of resonance. The current thesis focus on how this wave-function are modified. 
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Chapter 3 

The Feshbach resonance in two-body 
physics 

As discussed in Sec. 12.31 a two-particle interaction in a dilute system is often approximated by a pseudo- 
potential characterized with a s-wave scattering length . The drastic change of as obtained by tuning the 
energy difference between two channels through a magnetic field in the Feshbach resonance gives experimen- 
talists a rare ability to tune the interaction strength between two atoms. And this possibility is extremely 
useful to study BEC-BCS crossover where the interaction varies from weak to strong. 

Here we briefly review the Feshbach resonance in a two-body system. As discussed in Chapter [21 a 
hyperfine level is an eigenstate for an isolated single atom. However, when two atoms interact, most of the 
interaction comes from electrons with only negligible effects from nucleons . As a result, hyperfine levels 
are no longer true eigenstates of the two-body system. Nevertheless, hyperfine levels can serve as good 
approximated quantum numbers and we are going to call a pair of hyperfine indices a "channel" . Different 
channels in general have different interactions. They are decoupled at the lowest order. In a magnetic field, 
different channels differ in energy at threshold, where two atoms are infinitely away from each other, due to 
the Zeeman energy which is mostly determined by the electronic magnetic moment because the electronic 
magnetic moment is much larger than the nuclear magnetic moment. This energy difference is easy to tune 
through a magnetic field. 

When the mixing between channels are taken into consideration, the simple single-channel scattering 
becomes the multi-channel scattering. Especially, when the one channel's threshold is close to a bound-state 
in the other channel, the low-energy scattering property of that channel is dramatically altered. Its phase 
shift changes 27r; its s-wave scattering length blows to infinity and then jumps to the infinity of the 
opposite sign. This is essentially what happens in a Feshbach resonance, studied by Fano [Fano, 1961| and 
Feshbach [Feshbach, 1962| for nuclear and atomic physics in 1960s. Here we mostly follow the treatment by 
Leggett [Leggett, 2006[ (with some different symbols to comply with the rest of this thesis). Q 

^In order to be consistent with other parts of the thesis, we here use some different symbols comparing to the original works 
by Leggett [Leggett, 2006| . Here we list them, with the symbols from [Leggett, 2006] in parenthesis. U (= —V): open-channel 
interaction; V (= — Vc): closed-channel interaction; Y {=-g ■ /): inter-channel coupling; _E(, (= to): binding energy of the 
closed-channel bound state; Tc (= tq): potential range; r) {= + 5): the Zeeman energy difference between two channels; K 
(= k) see its definition in Ea l3.12l 
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The Hamiltonian for the coupled open- and cfosed- channel can be written as a 2 x 2 matrix 

/ 



H{t) 



\ 



-i^y^~U{r) -Y{r) 

-Y{r) -^v2+r,-y(r)^ 



(3.1) 



where the zero energy is taken as the energy of two atoms in the open-channel with infinite separation. 77 
is the difference in the Zeeman energies of the two channels. The first column (row) stands for the open- 
channel and the second column (row) stands for the closed-channel. All the interactions are short-range. 
For a s-wave solution, we have the radial part as 



(3.2) 




The coupled time-independent Schrodinger equations in the radial direction read as: 



2m 



-X" -Ux- Yxc = Ex (3.3) 



2mr 



Xc + VXc - Vxc 'Yx^ Exc (3.4) 



We now expand the closed-channel component Xc over the eigenstates of the isolated closed-channel 
Hamiltonian, Xc = "Ylii ^ifj^i, where (jji satisfies Schrodinger equation of the isolated closed-channel 



Vcf>, = -Ei;><j>i (3.5) 



We denote 0o the wave function in resonance and cq its coefRcient. Here we assume that the energy 
differences between eigenstates, 0i's, are larger than any other energy scales in the problem. This guarantees 
Co dominates any other q's. So, Xc ~ Co0o- Naively speaking, the resonance is expected to happen at the 
point where the closed-channel bound state has its energy exactly at the threshold of the open-channel. This 
leads us to introduce the relative detuning, 6 = 1] — Eh. By comparing Eq. 13.31 and Eq. 13. 5[ it is easy to 
show 

Xc = ^Jdr'roir')Y{r')x{r') (3.6) 
provided that 00 is normalized (for radial component). Inserting the expression of Xc into the Schrodinger 
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equation of the open-channel component Eq. 13.31 we get 

^2 „ ... 1 



X" -{U + E)x + ^ / K{r r')x{r')dr' = (3.7) 



2m,r E - S Jo 

where the kernel K{rr') is 

K{rr') EE (t>o{r)Ur')Y{r')Y{r) (3.8) 
The Schrodinger equation of the decoupled open-channel is 



2mr 



X" -iU + E)xa = (3.9) 



Note that Xoi''') and Yq(^) '^^^ A{1 — r/abg) , where abg is the background s-wave scattering length 



of the isolated open-channel 



Multiply Eq. 13.71 by xo and Eq. 13.91 by x, integrate from r = to a distance 



ro much larger than the potential range, Tc, subtract them, we find, 

dr[xo{r)x"{r) ~ Xo{'r)x{r)]^ ^ drxo{r)x{r) = J ~, dr dr' xo{r)K{rr')x{r') 

n Jo E — Jo Jo 

Using the Green's theorem on the first integral, we find 

2mrE r" 2mr/h^ 
Xo{ro)x' (ro) ~ x'o{ro)x{ro) + J, / drxo{r)x{r) ^ —^-^ dr dr'xo{r)K{rr')x{r') (3.10) 

« Jo E-S Jo Jo 

Here we can use the boundary condition by Xo(^) A{1 — r/a^g), x(^) A{1 — r/ag). Y{r) is a 

short-range interaction and thus K(r,r') only picks the short-range parts of x(?') and Xo(^)7 which yaries 
little for different detuning. Consequently, the R.H.S approaches a constant. For the scattering solution at 
E — 0, the last terml on the L.H.S. is zero, and we haye 

- - — = r dr r dr'xo{r)K{rr')x{r') (3.11) 

as abg d Jo Jo 

Contrary to the preyious guess, a^ does not diyerge at the point 6 = because of the original interaction in 
the open-channel, abg. We can introduce a quantity /C, the detuning where as diyerges, through the implicit 
equation (/C shows up in R.H.S as well) 



oo poo 



JC^-^^I dr I dr'xo{r)K{rr')x~s^^ir') (3.12) 



"'0 



■^Note that we are dealing with the internal wave function xo here, i.e. the wave-function in the region I, instead of the 
external wave function as in Sec. 12.31 therefore, the boundary condition for r — > in Sec. I2.3l actuallv corresponds the boundary 
condition r — > oo here. 
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And if we define tlie "real detuning" S — S — K., which is the detuning from the real resonant point, where 
as — > ±00, we have 

as{S) ^ abg (^1 + (3.13) 
This is consistent with the empirical formula of the Feshbach resonance 



as{B) = flbg 1 + 



AB 
B-B, 



(3.14) 



where Bq is the magnetic field at which the diverges, i.e., the resonant point. Comparing the above two 
equations, we see that AB = K.{d6 / dB)^^ . 

flj/abg 




-10 



-5 



0.5 

-0.5 
-1 
-1.5 



Figure 3.1: S-wave scattering length in Feshbach resonance 
The dashed line is j/ = abg. 
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(B-B„)/AB 



Let us now consider the bound state rj where E <Q. We define 



ab{E) 



(2m, 1^1)1/2 



(3.15) 



Here, we only study the bound state close to threshold with a binding energy much smaller than the binding 



^We wish to stress again the difference between a uncoupled closed-channel bound-state (l>i and a real bound state formed by 
atoms from both open- and closed-channel. The closed-channel bound state, <f>i, is the eigenstate of the isolated closed-channel 
Hamiltonian. It is not a real eigenstate for the full two-channel Hamiltonian. On the other hand, the full Hamiltonian has 
bound eigenstates (i.e. i? < 0) at large negative detuning (in the two-body context). The bound state solution for a two-channel 
eigenstate has components in both the open-channel and the closed-channel. The open-channel weight is often larger when close 
to the resonance; consequently, such bound-states are often called open-channel bound-states. Only at large negative detuning, 
the real two-channel bound state, mostly composed of close-channel component, coincides with the close-channel bound-state, 
(l>i to a large degree. 
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energy of the uncoupled closed-channel bound state, (jjQ, namely, \E\ <^ Eh. Therefore, we have Of, ^ Oc. 
Outside the range of potential rc, the wave function is proportional to e"''/"''. For Tc <C r ^ at,, this wave- 
function can be expanded as 1 — just as the Bethe-Peierls boundary condition (or the s-wave scattering 
wave function) we discussed in Chapter 12.31 Note that ab{E) is not identified as a priori. Through a 
procedure similar as previous, it is not hard to find 

flh + K, — E 

We have assumed the short-range part of the bound-state wave function x does not change much from the 
short-range part of the scattering state 3; therefore JC stays like a constant. Provided both S and \E\ are 
much smaller than /C, this yields 

^^ = ^^9J~~^ (3-17) 

It is not hard to see that indeed coincides with the s-wave scattering length as in Eq. I3.13l when \E\ <C 
Thus we will use ab and as interchangeably hereafter. This is actually an example of our discussion about 
the Bethe-Peierls boundary condition in Chapter 12.31 Using the concept of the Bethe-Peierls boundary 
condition, we do not need to distinguish the zero energy scattering state and the negative energy bound 
state, and we can arrive the same conclusion easily. 

Using Eqs. I3.15l and l3.171 it is easy to obtain an equation for E 

{\E\+5f -25c\E\=Q (3.18) 



where 5c is defined as 



h? /mralg 

And the solution for the negative detuning, i.e. 5 < 0, is 



(3.19) 



\E\=5c-5- V^c - (3.20) 



We can also calculate the "relative weight" on probability of the closed-channel for a normalized open-channel 

*This is guaranteed by the boundary condition of x{^) — > 1 — r/as, which fixes the normalization of the short-range 
part of the wave function (x?* — > 0) to be the same, namely 1. As a result, all the integral term over the short-range kernel 
fj drdr'A'(r, r')xo{^)x{^') gi™ the similar value as K at the resonant point, Eq. 13.121 



16 



component, Xn, of the wave function (/ Xn{r)'^dr = 1) 



(3.21) 



Comparing this with Eq. I3.12[ and assuming the short-range part of x„(r) does not differ from Xo(^) much 
except the normalization, wc find 



For \5\ ^ Sc IC, we have 



We can simphfy A further 



X-^^^^ (3.22) 



E-S^IC (3.23) 



|£;|Xi/2 



2mrabg {JCub) \26c) (3.24) 
This shows that Sc is a "characteristic" energy scale. When \S\ ^ 6c, Eq. 13.201 gives E 5 , and the above 



equation gives A — yj\5\/ {25c) 3> 1: this simply means that when the negative detuning is large, weight is 
predominently in the closed-channel; the real mixed bound state is a closed-channel bound state {4>o) slightly 
dressed with the open-channel component; therefore the binding energy is very close to the binding energy 
of the uncoupled closed-channel bound state, 0o. The often quoted relation between as and the binding 
energy, Eq. I3.15| is not valid here. On the contrary, when \5\ <^ Sc, Eq. 13.201 gives \E\ « S^/Sc <C \S\, 
and Eq. 13.241 then gives A « S/{^/25c) <C 1. Here, the closed-channel weight is much smaller than that 
of the open-channel; the real bound-state is essentially an open-channel affair with little dress-up from the 
closed-channel. Eq. 13.151 is valid in this case. 

When dealing with many-body problems, another important energy scale comes into play, namely, the 
Fermi energy, Ep. When Ep is much smaller than Sc, i.e. a broad resonance, the closed-channel has 
only negligible weight close to resonance, which is the situation of interest. In such a case, it is a good 
approximation to take the Feshbach resonance as only a knob to tweak the interaction in the open-channel. 
On the contrary, when Ep is close or even larger than Sc, i.e. a narrow resonance, the closed-channel weight 
can be significant clsoe to resonance; as a result, it is necessary to explicitly take the closed-channel into 
account in the many-body framework. 

One important remark about the above discussion is that whether a resonance is "broad" or "narrow" is 
a many-body concept and is relative. For one particular resonance with a fixed Sc, in a infinite dilute system 



^Note here the normaUzation on x{^) different from the most cases of this section, where x{^) normahzed by requiring 
X (1 — r/ai,) or e~^^'^i> here. The difference is an extra (a;,)"^/^ for the wave function and that contributes the extra 

factor in Eq. |3^ 
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(i.e. a two-body system) with a zero Fermi energy, Ep = <C (5c, it is a broad resonance. As the system 
becomes denser and denser, Ep increases and the resonance becomes less and less broad; finally, when Ep 
is in the order or even larger than Sc, the resonance becomes a narrow one. In the real experimental system 
of the ultracold alkali gases, the achievable density has a certain range, so does the Fermi energy range, Ep; 
therefore the "broadness" or "narrowness" of a resonance is not as flexible. 
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Chapter 4 

The single-channel BEC-BCS 
crossover 



In this chapter, we briefly review the BEC-BCS crossover in a single channel. The idea to describe the BEC 

and BCS on the same footing stems back several decades [Eagles, 1969|Leggett, 1980|Nozieres and Schmitt-Rink, 1985[ 

IRanderia, 1995] . The BCS theory can be understood as the following: fermions form "giant molecules" and 

those molecules then condense simultaneously. It is not hard to show that the BCS ansatz is equivalent to 

the coherent state of a two-body pair -0'!' — Cka|j.a|_]j. 

n("k + t^k4k4-k) |0) = ^exp Ck4k4-k) |0) (4-1) 

k k 

where Ck — (fk/wk). The size of these "giant molecules" is much larger than the interparticle distance. 
Moving away from the BCS end toward the BEC end, pairs shrink as the interparticle attraction becomes 
stronger. On the BEC side, two-fermion molecules are smaller than the interparticle distance and therefore 
a well-defined object. However, the binding energy now is much larger than the typical many-body energy 
scale and condensation does not happen at the same time when molecules form. Nevertheless, at low enough 
temperature, we can still consider the formation of molecules and their condensation at the same time, 
compatible with the the same framework used for BCS. 

In the single-channel BEC-BCS crossover model, one imagines a "magic" knob that can tune the inter- 
action strength along the crossover. The many-body fermion system sweeps from BCS of a fermionic atom 
system to BEC of diatomic molecules in response to the increase of attraction. This directly applies to 
the broad-resonance of the two-channel case as well, where the closed-channel weight is negligible and only 
serves to modify the effective interaction strength in the open-channel. We mostly follow the path-integral 
treatment by Randeria and the company Randeria, 1995pi!ngelbrecht et al., 1997[|Diener et al, 20081 . They 



have studied this problem with a path integral approach which is proved to be a nice tool for the problem 
due to its flexibility and readiness to be extended for the higher order fluctuations. In the next two chapters, 
this method will be adapted further for the two-channel model. 

We start with an attractive J-potcntial in the coordinate space. This potential is not equivalent to the 
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reduced pairing potential used in the original BCS work. The reduced pairing potential only couples particles 
of the opposite momentum and does not support simple form of Hubbard-Stratonovich transformation, which 
is essential to solve the problem in the path integral formulation. 

The Hamiltonian with the chemical potential of the system can be written as 

H-y.N^Y.j '^''^^^('•) ("i^' - - 9 J d'^r4(r)c|(r)c^(r)c^(r) (4.2) 

Introducing the quantum partition function 2 = J 'D{ip, ip) exp {—S[ip, ipf) , where tp) denotes the func- 
tional integral over all possible wave function ijj and ijj, and the action Slip^ip] can be written down from the 
Hamiltonian 



Slip^ip] = J dT I d^r 



1 



'ipair,T) [ dr - —V -/i] ipair,T) - ^^t (i", (r, t)'0^ (r, t) V'^ (r, t) 



(4.3) 



The fermion fields ipa and are two independent Grassmann variables. Notice that they are not complex 
conjugate to each other as in the usual operator language because complex conjugate is not a well-defined 
concept for Grassmann variables. 

This system can be solved with Hubbard-Stratonovich transformation. Introduce an auxiliary field 
(functional variable) A(r, r) coupled with a pair ^p■\{r, T)^pi{r, r). We write down first the Gaussian integral 
of A 

^ ' ' ■ - ^ (4.4) 



1 = J V{A,A) exp (^-^ J drd'^r AaJ 



Note that we absorb the extra constant of integration into the measure of I?(A, A). And with a shift of 
A(r, r) -> A(r,T) - gip^{r,T)i;i{r,T)), we havell 



exp ( g / drd'^rV^t^^V'^V't ) = / 2?(A, A) exp <{ - / drd'^r 



^AA - (AV^;Vt + AV^tV'-i) 



(4.5) 



Note that A(r,T) (or A(r, r)) comes from Grassmann fields ■0(r, r) (or ^(r,T)). Therefore, they are not 
related to each other as complex conjugate either. Nevertheless, at the mean field level or only at the phase 
fiuctuation around the mean field values, A and A are indeed complex conjugate. Consequently, we will just 
take A as normal bosonic field in the following and often simply treat A as A's complex conjugate. Now 



^ J 'D{A, A) ■ 1 is only a constant factor on partition function Z and has no effect on real physical quantity; therefore, we 
can take it as 1. (This is equivalent to divide the 2 by a constant) 
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the interaction term can be replaced. 



Z=J J 2?(A,A) 



exp ■ 



drd^r 



At the expense of introducing an auxihary field (A) which has contact-type coupling to the original field 
ip, we eliminate the four-field interaction term formally. A field is like a local potential for ■(/;, although this 
local potential has to be calculated from the original field self-consistently. Nevertheless, A couples to a 
pair of fermionic field ip, and thus it extracts a special degree of freedom from the field. When properly 
selected, this degree of freedom is highly non-trivial and has macroscopic importance, which serves as the 
"order parameter" for the system. The above formula for partition function is bilinear to ip, and we can 
rewrite it into a nicer form in Nambu spinor representation 



5- = 



(4.6) 



P(*,4') / 2?(A, A)exp 



drd'^r 



-AA- 



where 



V 




is known as the Gor'kov Green function. [Gq ] 



(p)i-i 



-Or 



2 m 



Ai, and [G^oY' 



-dr 



(4.7) 



(4.8) 



2m 



represent the non-interacting Green's functions of the particle and the hole respectively. 

Before going further, we would like to discuss one confusing point about the possible one-or-two indices 
for quantities such as G or A in Eq. 14.81 As a matrix, such a quantity has two indices (x, x') or {p,p'), which 
have no ambiguity in usage. On the other hand, there are often ambiguity when only one index x or p is 
used. In some cases, the one index means the relative value of the two indices. For example, an interaction, 
[/(ri,r2), normally only depends on the relative coordinate, r = ri — r2. So C/(r) means C/(ro + r, Tq). In 
other cases, especially common in the current thesis, the one index stands for its repetition. In this case, 
the difference is always zero. For example, the free Green's function, Go{p) stands for Go{p,p')S{p — p'). 
Similarly, the order parameter, A(a;), only couples to ip(x)ip{x) (Eq. 14. 5p . When used in the matrix context 
(Eq. 14. 8p . it means A{x)S{x — x'). Interestingly, its Fourier transformation in momentum space does not 
have the same properties. In fact, it means A{p,p') — A(p' —p). 
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Now action in Eq. 14.71 is bilinear to ^; so it can be integrated out formally and the partition function 
then only depends on the field A. 



And the action becomes 



^(A, A)exp 



5[A,A] 



drd'^r- Aa] -Indet^" 
9 J 



j drd'^r^AA^ -Indet^"^ 



(4.9) 



(4.10) 



Note that the determinant in Indet^^^ runs through both the normal coordinate space and 2x2 Nambu 
spinor space. The above formulas are exactly equivalent to the original partition function (action) in the 
fermion field (Eq. 14. 3p . It looks nice and compact. Nevertheless, Indet term is highly non-trivial and 
contains all the many-body physics. 



4.1 Mean field results 



The saddle point equation of Eq. (|4.9p gives the mean-field result of the system. First we need to find the 
derivative of Indet^"^. We notice the identity 



In det ^ = tr In A 



(4.11) 



and differential rule of the function "tr In" 



-^trln(r^) = tr(^-^^-i) 

)(pq 0(pq 



(4.12) 



Using the above relations, we can write the saddle equation of Eq. (|4.9p (differential with respect to A) as 



-A(r,T) - tr 



^(r,T,r,T) 



1 



0, 



(4.13) 



Here this matrix is in the Nambu spinor space. At the mean field level, we seek a tempo-spacial homogeneous 
solution of A(a;) = Aq. At this level, A{p) becomes a (5-function in the frequency-momentum space, and 
has non-zero elements only for two fermions with the same momentum. (Please See the discussion in the 
previous section about one vs. two indices. This is not generally true in other situations, as we show it 
when discussing collective modes in sec. 14.21) We can find the Gor'kov Green function from Eq. (|4.8|) in 
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momentum space at the mean-field level 

/ 



-Ao 




Go(p)W (4.14) 



Here p is the frequency- momentum, p = {ujn, p), and a;„ is the Matsubara frequency of Fermions. = eu — 



Ek — k /2r7i, Ep = + |Ao| . And the saddle point equation can be rewritten as 

-^° = ^E^fW (4.15) 

9 Vo ^ ujI + El 

Here T is the ternperature, and Vo is the volume in d-dimension. The summation of the Matsubara frequency 



can be evaluatec ^ 



and we find 



1 _ 1 V 1 - 2»/(igp) _ 1 V- tanh(igp/2r) 

g Vo ^ 2Ep 2Ep ^ ' ' 

where n / (e) is the fermi distribution function. This is exactly the famous gap equation obtained from other 
methods as well. On the other hand, Q^^ in Eq. (|4.8p is the inverse of the fermion-fermion correlation 
of 4*. In the mean field, Gq as Eq. (j4.14p can be diagonalized in the momentum space with a canonical 
(Bogoliubov) transformation. We can make an analytic continuation of zw„ — > w + 0+. Eq. (I4.14p then 
has poles {±Ep) where uj"^ — E^ = 0, which determine the spectrum of fermionic excitations. Indeed, in the 
BCS-likc states (/i > 0), the spectrum is gapped at A; while in the BEC-like states (/i < 0), the fermionic 
excitation starts from the molecule binding energy -^//i^ + A-^ « 

The summation in Eq. ()4.16p does not converges in 3D because the summand does not decreases fast 
enough. This is because our assumption of contact interaction breaks down for the scale smaller than real 
potential range Tc, i.e., the summation of momentum is capped at some high momentum A related to l/rc- 
Notice that in 3D, we have a similar relation that connect the bare potential 5 to a more physically observable 
quantity, the s-wave scattering length 



mVo 



^The summation of the Matsubara frequency of a function h{iuj„) is carried out by the normal trick. We multiply h{z) with 
the Fermi distribution function np{z), the summation is the sum of residuals at the imaginary axis of np{z). The contour 
can be deform into a contour over the rest of singular points of h{z). We just need to find the residuals of the total function 
h(z)np(z) over those singular points to find the Matsubara summation. However, due to zero temperature, the np{z) is only 
nonzero at the negative singular points of h(z), — -Ek in this case. (The other singular point, Ek; gives np{Ei^) = for zero 
temperature.) 
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Here Vo is the total volume. We can renormalize Eq. 14.161 with this relation 



mVo 
47ras 



E 



iimh{Ek/2T) 1 
2Ek 2ek 



(4.18) 



Now the gap equation has proper decay in high momentum and no artificial cutoff is necessary. There are 
two unknown parameters, and A, in the equation. We need another equation in order to pin them down. 
To complement the gap equation, we can introduce the number equation, N = —dfl/dfi. At the saddle 
point, the thermodynamic potential is flo = ^[Ao]//?, and we have the number equation 



Similarly the summation (due to the trace) over the Mastubara frequency can be evaluated and we have the 
number equation 



This equation has no divergence at high momentum. The number equation and the renormalized gap 
equation Eq. (I4.18P compose the implicit equations for two unknown parameters, gap A and chemical 
potential fi. It is not hard to find the zero temperature analytic result at both ends. At the BCS end 
(l/kpUs — )■ — oo), we obtain jj, w Ep and A oc exp(— 7r/2fcF \o-s\)', at the BEC end {l/kpCLs ^ +oo), /i = 



significance. In the more general crossover region, these two equations can only be solved numerically. 
They have no singularity in the whole region, which indicates it is a crossover instead of any simple phase 
transition. Please see Fig. l4.1l for detail. 

4.2 Gaussian fluctuation and collective modes 

We can expand the partition function Eq. (14.91) around the mean-field value, A(r,T) — Aq + 6'(r, r). The 
linear order of the expansion is zero because Aq is the saddle point. The next order gives us the bilinear 
terms on 0, i.e., correlation of bosonic fields A. Note that here the Hamiltonian only has a contact- type 
potential, therefore it cannot cover the situation of a charged system where long-range Columnb interaction 
cannot be neglected. We limit ourselves to the neutual case. Nevertheless, it is conceivable that a more 
realistic short-range potential only renormalizes some parameters in the following calculation while leaves 
the qualitative result unmodified. 





(4.19) 



fi^ /2mal, i.e. half of the binding energy of a molecule, while A cx n^^^ag 



no longer has much physical 
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Figure 4.1: The chemical potential ^ and gap A in the mean field level over crossover 
All quantities in the unit of energy {fj,, A) are rescaled with the Fermi energy Ef and the s-wave scattering length 

as is rescaled with l/kp- 



Notice that we can expand the second term in Eq. (14.91) for G ^ — Gq^ + K 

trlnG"^ = trlnGo"Vtr(Goi^) - itr(Go^Go^) + ••• (4.20) 

In our case, 




(4.21) 



Here the linear terms oi K or 9 {6*) are zero as the saddle point condition. To the second order, the action 
is 

S[Ao, 9, 9*] = S[Ao] + ^tY(KK) + UiiG^kG^k) (4.22) 
2.9 2 

Write the last term into the momentum representation 

tr(Go-ft:GoX) = 5Z (Go(p)/^,Go(p - q)K_q) (4.23) 

Notice that the second "Tr" and following "Tr" in this section only runs in Nambu spinor space and q = 
{(\,qi), p — (p,Pn) are all four momentum, where qi is the bosonic Matsubara frequency while p„ is the 
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fermionic Matsubara frequency. 



e„ 



And we remember that Go{p) = Gop,p If we introduce a new vector 



e{q) 



0^1) = [e* 0^ 



the action can be rewritten into a more compact form 



S[Ao,e,e*] = S[Ao] + [^t(q)M(g)(^(g)] 



(4.24) 



(4.25) 



(4.26) 



Notice that we can always choose a real Aq and therefore Gq 12 (p) = Gq 21 (p), we have 



M,,, = M(g) = 



l + EpGo 11 (p) Go 22 (p- 9) EpGo i2(p)Go i2ip-q) 



\ 



(4.27) 



Ep Go i2{p)Go i2{p -q) ^ + Ep Go nip - q)Go 22{p)j 
The summation over the (fermionic) Matsubara frequency of Pn can be carried out at zero temperature 



Mii(g) = M22{-q) 

= - + ^ Go ii(p)Go 22{p - q) 

^ P,Pn 



9 ^\iqi-E-E' iqi+E + E' 



Mi2{q) = M2i{q) 

= ^Go i2(p)Go 12 (p - q) 



P,Pn 



= > uvu V — —7 — — 7 

^ \iqi+E + E' iqi-E-E' 



(4.28) 



(4.29) 



where u = Up, v = Vp, E = Ep and u' = Up-q, v' = fk-q, E' = E'k-q- Wk, ^^k, E^ are defined as usual BCS 
literature. 

^^ = l-«^ = ^(l-;|) (4.30) 
The G^^^ = M"^ is the correlation function of 9 (or A) and its poles give the spectrum of collective modes 
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as every 9q (or Ag) involves many fermions moving in a coherent manner. So the spectrum of cohective 
modes can be determined by finding poles of G'*^-*, det Af(a;,q) = 0, after we analytically continue for the 
frequency iqi ^ lo + iO+. 

For low energy modes, where w, |q|^ both are much smaller than min{i?k} = Aq (or -^//i^ + A^ for 
/i < 0), we can expand M with lo and q. The lowest order has the form uj Ri cq, which suggests a sound 
wave as expected for any Goldstone mode. At BCS side, c = vp/\/3, where vp is the Fermi velocity. 
This coincides with the famous Anderson-Bogoliubov mode. At the BEG side, we get = A^/8m|/i| = 
v'jp{kpas) I'i-K = ATrnBCLB /iTiB, which fits the low momentum part of Bogoliubov spectrum of bosons gas. Here 
ms — 'im is the molecule mass, ub = n/2 is the molecule density and as = is the inferred interaction 
between molecules. This value differs from the result of more accurate calculation from the few-body theory, 
as = Q.Qas [Petrov et al., 2004] , which indicates the possible deficiency of the current theory. 



4.3 An alternative method to invert the Green's function 



In the above section, we inverted the Gor'kov green function matrix Eqs. (j4.81 I4.14[) directly and it is 
not hard to do as a 2 x 2 matrix in the momentum space. Alternatively, we can use a different approach 
which proves to be more convenient in the two-channel problem. First, we diagonalize Q^^ with a unitary 
transformation T, in the momentum space 



ii^n - £,k A 

A iujn + Cfc 



= T^BT 



(4.31) 



It is easy to show that such T and B satisfying above equation are 



T = 


1 


Vk 




iuJn + Ek 






Uk. 




^ iuJn - 



(4.32) 



Eh 



where u^(w^) — ^(1 ± ^k/Ek) and E^ — ^/^+A^ are conventionally defined quantities in the BGS theory. 
Actually, this transformation is nothing but the Bogoliubov canonical transformation, and the B matrix 
simply describes the spectrum of the fermionic quasi-particles. Now it is easy to invert Q^^ 



g = T^'B^^T 



(4.33) 
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Green's function Q takes a more conventional form A/ (icUniEk) here without any dependency on frequency in 
nominator as Eq. (j4.14p . Matsubara frequency summation over Go(fc) in the mean-field and Go{k)Go{k + q) 
in the Gaussian order are then easier to perform as in text-book. 
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Chapter 5 

The two-channel three-species 
many-body model, the mean field 

For the narrow Feshbach resonance, atoms have considerable weight in the closed-channel and the Pauli 
exclusion between two channels cannot be neglected. A many-body framework needs to include both chan- 
nels. Before diving into the detailed calculation, let us make some rough estimates about scales in order 
to build some intuition. The same problem at the two-body level is well understood as briefed in Chapter 
[21 We estimate how much each channel would change in a many-body system. Compared to a two-body 
system, the two-body correlation of a (cold) many-body fermion system in the momentum representation is 
modified mostly in low momentum, i.e. around or below the Fermi momentum, while staying just like the 
two-body wave-function intact in high momentum. The open-channel component of the two-body scattering 
wave-function is like a zero-energy free wave (fc = 0) plus a small short-range kernel of size r^- In the 
momentum space, this wave function is like a J-function at fc = with a small tail in high momentum. 
The occupation number in the first available level, fc = 0, is close to one (for two hyperfine spins) and 
this level cannot accommodate one more pair. The next pair has to occupy the next available level in k 
instead of the lowest energy two-body level, fc = 0. This suggests that the open-channel requires a full 
many-body treatment. The situation is quite different in the closed-channel. One crucial assumption is that 
the closed-channel bound state is much smaller than the interparticle distance in the real representation. A 
closed-channel bound state, 0Oi in the momentum, spreads its most weight in the range [0, l/cc]. The weight 
in each momentum level is so tiny that it is even smaller than 1 when it is multiplied by total atoms number 
N , N \(j)k\^ <C 1. It only has a very small fraction of its weight within the Fermi energy range (< kp)- Hence, 
the closed-channel only has a small overlap with atoms in the open-channel, as well as with other atoms 
in the closed-channel. The smallness of this overlap ensures that many-body effects in the closed-channel 
can be treated perturbatively. As discussed in Sec. 12.31 and Appendix IB. 31 the high momentum part of the 
two-body correlation just follows the two-body wave function; therefore, we can write two-body correlation 
as 

/ik ~ 0ok/(k) (5.1) 
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with /(k) = 1 for kp. /(k) deviates from 1 in low energy and represents all the many-body correction. 

In addition to the two channel description, the model needs to be explicitly expressed with the three 
hyperfine species at least in some part in order to address the inter-channel Pauli exclusion between two 
channels due to the common species. This effect has received little attention in theoretical research and a 
model as such is a useful addition to our knowledge of the BEC-BCS crossover and the Feshbach resonance. 
Nevertheless, as we just discussed, and will illustrate more quantitatively later, the effects of the inter- 
channel Pauli exclusion between the two channels is relatively minor and can be treated perturbatively. 
Consequently, it is not necessarily to carry the three species description all through the calculation. The 
description of two channels and the description of three species are used in different phases of the calculation 
to solve the problem. 



5.1 The extremely narrow resonance 

Before the quantitative model, we discuss qualitatively a simple yet revealing case, the extremely narrow 
resonance, where the inter-channel coupling approaches zero, 1" — !■ 0. In this case, two channels are almost 
independent to each other except sharing the same chemical potentials. From the two-body discussion in 
Chapter [3] (Eqs. 13.81 13.121 I3.19p . the characteristic width of the resonance 5c is zero in this case, which 
indicates that the resonance is always narrow no matter how dilute the system is. 
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Figure 5.1: Extremely narrow resonance 
The shaded area is occupied by atoms. 
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Here we assume that the atoms in the open-channel maintains a relatively clear Fermi surface, like the 
Fermi liquid or the BCS. The many body situation is quite straightforward. In general, three different 
situations exist, as Figs. IS.la) ISTbl and lSTcl When the closed-channel bound state level is above the Fermi 
sea {Ep < 6, Fig. 15. lap , all atoms are in the open-channel and no atoms are in the closed-channel. Therefore, 
it is a genuine single-channel problem (only open-channel) with the original bare open-channel interaction. 
There is no need to considering the inter-channel Pauli exclusion. When the closed-channel bound state 
level is below the Fermi sea {6 < 0, Fig. IS.lcp . all atoms are in the closed-channel and no atoms are in 
the open-channel. It is a simple di-atom molecule gas. There is no need to consider the inter-channel Pauli 



30 



exclusion either. The interesting situation is when the closed-channel bound state level is in the Fermi sea 
{0 < S < Ef, Fig. [5lb)) . The Fermi sea is then filled from the bottom (zero energy) all the way to the 
closed-channel bound state level in the open-channel. Then all the rest atoms go into the closed-channel. 
We do need to consider the inter-channel Pauli exclusion because atoms exist in both channel. 

We can imagine the two channels are decoupled and each has iV*^") {N^'^'> atoms. It is not hard to 
find the wave function in each case. The inter-channel Pauli-exclusion effect should be in the order of the 
overlap of these two wave functions. In the open-channel, the momentum space is filled with occupation 
1 from the bottom up to the chemical potential fi (or 5). In the closed-channel, as we discussed in the 
introduction, the occupation is roughly 1/Eb at the low momentum (< Ep), where Eh is the binding energy 
of the closed-channel bound state. So this effect should be roughly equal Ep/Ef,. 



5.2 Model set-up and the Hubbard- St rat onovich transformation 

One way to approach this problem is to extend BCS ansatz and then find the set of parameters that optimizes 
the free energy. We can write down a three species BCS ansatz as 

I*) = n ("k + v^albl^ + Wi^ald^) |0) (5.2) 

k 

with normalization, |uk|^ + l^^kl^ + — 1. Here "a" is the common species of two channels, (a, &) is 

the open channel and (a,c) is the closed channel. |0) is the vacuum state of fermions. This ansatz can be 
derived from a many-body coherent state of the two-body state as we did in the single- channel (Eq. 14.11) . 



I^f) oc exp 



|0) (5.3) 



0(^6) and 0^°'^) are related to Uk, I'k and as 

^i^"^ = ^ (5.4) 
= ^ (5.5) 

We can then proceed to find the set of (uk, ^k, w^) that optimizes the free energy. This method offers good 
intuition of the state with an explicit connection to two-body wave functions, ^f"'') and . However, it is 

not easy to find these parameters from optimization process. Furthermore, it is hard to extend this method 
beyond the mean-field to study phenomena such as collective modes. We briefly discuss this approach in 
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AppendixlX] Instead, in this chapter, we use the path integral method, which turns out to be more convenient 
for tlris problem. The Hubbard-Stratonovich transformation provides a powerful tool for studying non-trivial 
degrees of freedom (order parameters) in the system. It is more or less equivalent to other approaches at 
the mean-filed level. But it has great advantage to be easily extended to explore the fluctuation over the 
mean- field result. 

lian ; 



For a two-channel problem, we write down the Hamiltonian [U as 

1 



[Y-il:air)^pb{r)^cir)^a{r) + h.c. 



(5.6) 



Here r]j is the Zeeman energy of the specific hyperfine species, a, b, c stands for the three hyperfine species 
defined as same as in the ansatz approach (Eq. 15.21) . All the interactions {U, V, Y), are contact type, this 
simplifies the Hubbard-Stratonovich transformation considerably. It is plausible as we only study low-energy 
phenomena for the short-range potential. 

A contact interaction is fiat in momentum space. In other words, it pairs not only zero center-of-mass 
momentum pairs, but also finite center-of-mass momentum pairs. This is different from the paring interaction 



used in the original BCS work Bardeen et al., 1957 , which only pairs zero center-of-mass momentum, i.e. 
opposite-momentum, pairs. Nevertheless, the saddle point (mean-field) solution settles at the zero center-of- 
mass momentum pairing; therefore the mean-field solution coincides with that of variation method derived 
from pairing only opposite momentum atoms. On the other hand, the collective mode (fluctuation of order 
parameters) emerges naturally from the included non-zero center-of-mass momentum interaction. We do 
not need to introduce more general interaction terms beyond the simple opposite momentum pairing (e.g. 
the Coulomb interaction in [Anderson, 1958| ) to study collective modes. 

Three different types of Hilbert spaces are used in this chapter. The flrst one is the (infinite dimension) 
coordinate space and its reciprocal momentum space. The second one is the (3 dimension) space of hyperfine 
spin, a, b, c. These two are atom-based. The full many-body Hilbert space is iV-power (also direct-product 
type) of the direct product of these two. As discussed previously, in some cases, instead of a 3-dimension 
hyperfine spin space, we use the third type of spaces, open- and closed-channel, (a, 6), (a,c) (2 dimension). 
Strictly speaking, this one is just a subspace of the direct product of a pair hyperfine spin spaces. Neverthe- 
less, it is sufficient within our model because only these two combinations (channels) of a pair are considered. 



^Here and hereafter in the chapter, we take h = 1. 
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In principle, fermion fields ip and ip (and A and A defined according tojp and ip), botli Grassmann muiibers, 
are independent to each other and not related as complex conjugato This is marked by using (bar) 
instead of the normally used "^"(dagger) sign. In this chapter, a 3 x 3 matrix always refer to hyperfine spins; 
while a 2 X 2 matrix always refers to open- and closed-channel. And f sign is reserved only for hermitian 
conjugate of these two types of matrices. 

Introduce two channels into the vector form, {ijjip) is a column vector and {ijjip) is a row vector. 



(VV) 




The two-body interaction can then be written as a (2 x 2) hermitian matrix U in the channel space 



U = 



U Y 



Y* V 



(5.7) 



We can now write the Hamiltonian in a more compact form 



j—a,b.c 



2m 



The finite-temperature action is 



2m 



(5. 



(5.9) 



Similar as in Chapter 31 we can perform the Hubbard-Stratonovich transformation here. Introduce aux- 
iliary fields (functional variables), (Ai, A2), as a 2-component vector and start from the "fat identity", where 
all the integral constant is absorbed into the measure of functional integral of I?(A, A) [Altland and Simons, 2010] . 




(5.10) 



At = (Ai,A2) A = 
here x is four-coordinate, (r,T), J dx = Jq dr J d'^r. 

■^Complex conjugate is not a well-defined concept for Grassmann algebra 
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We can make a shift in A 



(5.11) 



Write it explicitly into the matrix form 

/ 




Ai, A2 



Ai,A2 



U Y* 
Y V 



First of all, notice that the new auxiliary field A is related to the mixture of two channels by a 2 x 2 matrix 
U. Also note that A.^ is not the complex conjugate of A^ in general as they are related to Grassmann 
fields 'ip'4' ^-iid V'V^ which are not complex conjugate to each other. But it can be verified later that, to the 
mean-field (saddle point) level as well as to the simple (phase-fluctuation) Gaussian level about collective 
modes, A^ and A^ are indeed complex conjugate (and real for saddle point). Hence, for the simplicity, we 
will treat them as such hereinafter and use A and A* more or less arbitrarily. The other point to notice is 
that, A(r, r) carries the the same coordinates as ip{r,T)ip{r,T) (or frequency-momentum coordinates in the 
reciprocal space) because there is only the contact interaction. Now the "fat identity" Eq. 15. 101 becomes 



1 = / V{Aj,Aj) exp { - / dx[A''U-^A - (V'Vi)A - A(V'V') + (V'V')t/(V'V')]} 



(5.12) 



The above equation use the fact U is hermitian, so WU ^ — U — I. It can be rearranged 



as 



dx{ip'il;)U{ip'iij)] = j V{A, A) exp{- J dx[A^U-^A - {4>^p)A - A(V'V')]} (5.13) 



exp I 



This is now ready to be applied to the original action in Eq. [57 



+ rij)^j + [A^tj-^A - {i)i>)A - A(V'V')] 



(5.14) 
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We can introduce a spinor similar to the Nambu spinor representation in the single-channel superconductivity. 



* = 



(5.15) 



The action can then be rewritten in a more compact form with respect to ^ and 4" 



(5.16) 



where 



V 



Ai 

A2 



Ai 



A2 



1 Y72 



(5.17) 



(5.18) 



2^V2-M + ^cy 

Rewriting in the frequency- momentum space, we find Q^^ decoupled in frequency and momentum. 

UuJn -ik-ria Ai A2 ^ 

Ai ILOn + ik+ rib 

^ A2 iuJn + ^fe + Vc J 

herqj = -^k^ ^ A*- The diagonal elements of the second column/row (6) and the third column/row (c) 
corresponding to negative energy, because of the particular opposite choice in the spinor {^2'^2 = '^b'^b and 
^3^3 = The non-diagonal elements mix V'a with V'b,c and therefore lead to a number- non-conserved 

theory. 

If we in addition assume rja — rjb ~ 0, rjc — t] (i.e., use rj as the absolute Zeeman energy difference of two 
channels) , in frequency-momentum space. 



(5.19) 







-Ai -A2 ^ 


g ^ = iujni - 


-Ai 


-6 




^-A2 


-{ik+v)j 



^In principle, there are two chemical potentials, fia and /ij, ^ because a does not convert to b or c. But we omit this for 
simplicity and absorb all the difference into 77;. 
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The action in Eq. 15. 161 is bilinear to 'I' and we can integrate out 'I' (^') formally 

S{A,A) = j dx[hU-^A-iT\ng-^'^ (5.20) 

Note that at this stage A(r, r) is not necessarily homogeneous in space or pseudo-time as in the mean- field 
result. 

5.3 Diagonalization of the Green's function 

Eq. (|5.20l) looks fairly simple and compact. Nevertheless, it has all the physics in it and is not as simple 
as it looks. The major problem comes from the term trln^^^, which includes logarithm and trace over an 
infinite-dimension matrix. All these operations are fairly straight-forward if we can diagonalize the Green's 
function (or its inverse) in a proper basis. It is not hard to see Q^^ is already decoupled in the frequency- 
momentum space (Eqs. I5.18[ [5T9)) . It is however mixed in the 3x3 hyperfine-spin space. The rest of this 
section is dedicated to diagonalize this 3x3 matrix in the hyperfine-species space. Note that in principle, the 
following discussion is not limited for constant A, but also applies to inhomogeneous A(a;„, k) as well because 
the 3x3 hyperfine space is independent to the frequency-momentum space. Here we use the approach in 
Sec. 14.31 to diagonalize it. In current problem, we need to diagonalize a 3 x 3 matrix (Eq. 15.191) . in other 
words, we need to figure out the Bogoliubov canonical transformation over which the Hamiltonian/action 
is diagonalized. Eigen-problem of the 3x3 matrix involves solving a cubic equation. An exact solution 
exists in principle. However, it offers little intuition to write down the exact result. Instead, the spectrum 
from the broad-resonance, where the only effect of the closed-channel is to modify the effective interaction 
of the open-channel, serves a reasonable lowest order approximation. We proceed to find the next order of 
correction over it (See Appendix IB. II and IB. 51) . 

We can break down the unitary transformation into two steps T and L. 

i3.„,k = LlTlGzl^nik (5.21) 

Here is the diagonal matrix; T and L are both unitary transformation. We take T as the canonical 
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transformation at the broad resonance, i.e., when we can ignore the inter-channel Pauli exclusion. 











-Vk 


Uk 




[ 






where Uk and Vk are defined in a similar fashion as in the single-channel BCS problem 



(5.22) 



2 — 1 2 _ 1 

= 1 - 1iv = - I 1 - — 



(5.23) 
(5.24) 



Note that here does not carry the physical meaning of the occupation number of the (open-channel) 
atoms, and i?k does not stand for fermionic excitation spectrum as in Chapter 2] or Appendix [X] They carry 
such meaning only at the broad resonance. In the narrow resonance, as we currently discuss, they are the 
zeroth order approximates of such quantities. 

In the broad resonance, the closed-channel can be integrated out at the two-body level and only the BCS 
pairing in the open channel needs to be considered at the many-body level. Matrix T, however, is enough 
to diagonalize and L is simply an identity matrix. In the narrow resonance, T cannot diagonalize 
because of the inter-channel Pauli exclusion between channels. Consequently, L stands the extra correction 
due to Pauli exclusion in the canonical transformation. Apply T onto G~^, we have 



/ 



TlGZ'MTk=z^nI- 



\Uk 



Ek 

+Ek 

A2 WfcA2 



MfcA2 ^ 
WfeA2 



(5.25) 



We regard the off-diagonal elements as perturbation because we only seek the solution around the BCS wave 
function (T transform). Introduce a dimensionless scale C, 



c 



A2 
^2 

Airy 



(5.26) 



Here both Ai and A2 are their mean-field (saddle point) values. It can be verified that C ^ 1 (See Appendix 
IB.5p . This matrix can then be diagonalized with the unitary transformation Lk within the first order of C, 
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(see Appendix IB.ll for details of calculation. 



Elk 






















-E, 



3k 



(5.27) 



The dispersion spectrum of fermions is 



E- 



Ik 



E' 



2k 



Eu 



7ik w -Ek + 



■ 72k ~ Eu 



^2"k 
^2^k 



E3W = ^k + ?7 + 73k ~ ?k + - 



£^k + UkC 



Sk - «kC 



Ck + j/ 



Ai^Ei^ + ulAiC 



V 



^2 



2(a + V) 



ek + V 



Ai « Sk - f^AiC 

2 Ck + 7? 2 



(5.28) 
(5.29) 
(5.30) 



The last step of each equations is valid only for low momentum {'^ kp) where 77 ^ ^k- Here we choose the 
sign convention to make -©1^2.3 positive in their zeroth order. Similarly as in Eq. 15.181 the second and third 
diagonal elements are negative because in the spinor representation, we choose ipb and ipc for ^E", which gives 
extra negative signs for quantity such as 5'^' in Eq. [51JI These negative signs disappear when we restore 
them to the normal order ^■0 in Sec. 16.11 

One interesting feature of this solution is that the corrections do not disappear for zero inter-channel 
coupling, Y — 0, which we discussed qualitatively in Sec. 15.11 This is because the corrections are due 
to the inter-channel Pauli exclusion between two channels, which continues to exist even when there is no 
inter-channel coupling. In fact, the inter-channel coupling, Y, contributing mostly energetically, modifies 
the effective interaction in the open-channel and affect the open-channel order parameters, Ai, greatly. By 
taking the broad resonance result as the zeroth order, we have taken effects of Y into consideration. On 
the other hand, the inter-channel Pauli exclusion or statistics is left out as the higher order correction as 
illustrated above. Furthermore, from Appendix IB. 51 (Eq. IB.28p . we know C ~ it is just the square root 
of our estimate Ep/Eb {Eb — h^K^ /2m) in Sec. 15.11 

The extra factor of unitary transformation is 





A1A2 



-JpT- Uk 



v 





-Wk 



Vk 





A2 



I + 6k 



4 



/-4 



(5.31) 



^Recall that jpi and tpi are Grassmann fields; thereore 4'i'4'i = —^i'4'i- 
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Here we use Uk'i'k — Ai/2ii^k- Note that L and L''' are unitary only to the first order of Ai/rj. Now it is easy 
to express the Green's function as 

G.„,k = T^L^B-^^lItI (5.32) 

This is ready to be expanded over the perturbation in order of ( or Ai/ij. It is easy to see that all uin 
dependence concentrates on -Bi^„,k, which is linear in ujn and simplifies the Matsubara frequency summation 
considerably. 

G.„,k « T^BZ'^^^Tl + T^5^Bzl^Tl - T^BZlJ^.Tl ^ G^^l^ + G^^^ (5.33a) 
G^'Im^T^BzI^tI (5.33b) 
Gi'ik = TMZl^Tl - T^BZlJ^Tl (5.33c) 



5.4 Mean field equations 



Use the same techniques for derivatives as Eq. (|4.12l) . we derive two saddle point equations for Ai and A2 
from Eq. ((5?20)) . 



5Ai 
8 



(;7-i)iiAi + (;7-i)2iA2-tr 
(;7-i)i2Ai + ([/-i)22A2-tr 



Go • 


/o 1 o\ 


000 


Vo 0/. 


Go • 


/o iv 


000 


Vo 0/ 



(5.34) 
(5.35) 



Taking A as real and constant ^ we can find the mean field result. Eq. (|5.19p can be inverted to get 
G. The inversion is quite tedious, but fortunately, we only need two elements of the G matrix (Gg (21) and 
Go (31)). The final mean-field equations are (For simplicity, both A^'s are taken as real. Please see Appendix 
R2]for detail) 

/ai\ ( JJ y\ ( /7i,.\ 

(5.36) 



where 





^Ik = (V-'a.-kV'fc.+k) = Ai 
/l2k = (^a.-kV'c.+k) = A2 



£^ik + a + »? 



{El k + k) (£"1 k + ^'s k) 
^Ik + Ck 

(£;ik + -B2k)(i;ik + -B3k) 



(5.37) 
(5.38) 



^When Y is not real, Ai and A2 cannot be both real even at the mean field level. Nevertheless, we can require one real, 
then the other will have a phase just to compensate the phase in Y . The final conclusion can be verified still valid. 
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Comparing Eq. 15.361 with the gap equation for the single-channel crossover, we see that Ai^2 are di- 
rect counterpart of the order parameter A in the single-channel problem, /iik and /i2k are the equal- 
time expectation of the anomalous Green's function for the open- and closed-channel. On the other hand, 
they correspond the macroscopic eigen-function of the two-body density matrix as described by Zhang and 
Leggett [Zhang, 2009[|Zhang and Leggett, 2009| (see Sec. 12. 3p . For the purpose of many calculations, they 
are just like the two-body wave function. And indeed, they coincide with the two-body wave function at 
high momentum. We will reference them often as the "two-body correlation" in the following. If we simply 
take Ei to the lowest order of and ignore the closed-channel, it is easy to identify /iik ~ Ai/(2£'k) as the 
many-body wave function in the single channel BEC-BCS crossover problem. 

At high-momentum, both hn^ and ft,2k behave as 1/ek which makes the summation (or the converted 
integral) diverges in 3D. This divergence can be mitigated by setting a high-momentum cutoff in integral or 
recognizing the decay of interaction in high momentum. It is important to recognize that this divergence 
is not unique in many-body and exists in the two-body physics, where high-momentum component of wave 
function is 1 /ek asymptotically as well. In the next section, we proceed to remove divergence of summation 
in ft-ik and ft,2k by noting the same divergence in the two-body wave function and manipulating accordingly. 

5.4.1 Renormalization of the mean field equations 

The mean-field equations (Eq. I5.36P can be rewritten as 

^ip = Yl f^pk/Mk + ^Pk'iSk (5.39) 

k k 

A2p = J2 ^Pk'^lk + J2 ^Pk'^Sk (5.40) 
k k 

The first thing to notice is that we restore the momentum dependence of the interaction as well as of Aip 
and A2p. Here Aip and A2p vary slowly in low momentum and it is a good approximation to take them 
as constant if wc only study the low momentum properties. We introduce the momentum-dependence in 
order to gain proper convergence in summation in high momentum. As pointed out previously, both hn^ 
and /i2k approach l/ck in high momentum asymptotically. This makes all summations diverge at high 
momentum when they are converted to integral in 3D if we take the interaction as contact and pull them 
out of the summation. However, we note that the real potential is not the contact-type and decays in high 
momentum. One way to remove the divergence is to introduce a cutoff at high-momentum while keeping 
the interaction coefficient constant. However, an arbitrary cutoff is undesirable for a theory. Alternatively, 
we can remove the divergence by restoring the momentum dependence in the interaction and recognizing 
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its decay at high-momentum. The problem here is that those microscopic bare interactions are hard to 
pin down. More readily available and observable is the two-body low-energy effective scattering matrix, 
T, which effectively integrates out the high-momentum component of the bare interaction. In the single 
channel problem, (Sec. 14. ip . we introduce a relation involving the two-body s-wave scattering length Os, (or 
To), Eq. (j4.17p . which has the similar divergence in high momentum. We subtract it from the summation 
and the remaining difference does not diverge in high momentum. A cutoff or bare momentum-dependent 
interaction is then no longer necessary. In principle, we can do the same in two channels and introduce the 
two-channel equivalent of the s-wave scattering length, a 2 x 2 scattering amplitude matrix. However, this 
2x2 matrix is awkward in definition and not clear for the physical meaning. Therefore, we adopt a two-step 
process with better physical intuition behind it. 

Besides the assumption of short-range potentials, we make another assumption: the two-body closed- 
channel bound-state at resonance, (j)o, is much smaller in size (oc) than the interparticle distance (ao), 
although it is larger than the potential range (re), Vc ac ^ ag. The problem would be a genuine 3-species 
problem and requires other techniques if the closed-channel bound-state is about the size or even larger than 
the interparticle distance. As discussed at the beginning of this chapter, this assumption guarantees that the 
closed-channel two-body correlation /i2k follows its two-body bound-state wave function in high momentum, 
and we can write it as 

hk 00k/ (k) dUD 

with the factor /(k) encapsulating all the many-body effects. This replacement also helps to cure the 
divergence involving /i2k because 0ok decays fast enough in high momentum and gives no singularity when 
integrated with the momentum-dependent interaction. After /i2k is taken care of, we proceed to remove 
the singularity of integral involving the open-channel two-body correlation /iik using the regular method by 
subtracting the open-channel low-energy two-body interacting kernel, or the open-channel s-wave scattering 
length, ai°\ which involves the same high- momentum divergence as ft-ik. 

Let us put the above qualitative statement into concrete mathematical formulas. As discussed before, we 
start with the closed-channel correlation ft,2k (I5q- I5.38p . In the lowest order of we have Ei^^^E2k ^ ^k' 
/i2k can be rewritten into such form 

^^^^ (Eik+Esk) ^ 2gk+r) j'^^^ same as two-body wave function (po besides normalization (when 
k <C k); while the extra factor is the many-body factor /(k) in Eq. I5.1l which describes the Pauli exclusion 
between two channels. At the low-momentum, the open-channel weight is large, the phase space left for 
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the closed-channel is limited, which is shown mathematically as a small factor. At higher momentum 
{kp <C fc ^ l/^'c), there is almost nothing in the open-channel, indicated by ~ 1 and thus the closed- 
channel wave-function just follows its two-body counterpart with a different normalization factor. It is not 
hard to see that A2 is closely related to the ^Hntegrated contact intensity^\ C, in Tan's work about univer- 
sality ( [Tan, 2008a[|Tan, 2008b| ). In his work. Tan concluded that the high-end of the relative-momentum 
distribution asymptotically approaches C /k^. Note that the high-end in his paper means momentum lower 
than 1/rc, but higher than any other scales (l/a^, 1/ao). In such scale, « 1, C = for the closed- 
channel. At even higher momentum (> l/^c), the mean- field / saddle-point solution no longer applies. And 
the two-body correlation, /12k, should just follow two-body wave function (See Sec. 12. 3p . By all the above 
arguments, we replace (^Ei^e^i^) '^ith ^ normalization factor, a, and the two-body wave function 0ok in Eq. 

ESI 

^2k = a</>ok'"k (5-42) 
where is the normalized solution of two-body Schrodinger equation, |</'ok|^ = 1- 

- £;^°Vop = 2ep(/)op - ^pk-^ok (5.43) 

k 

By replacing /12k with the two-body wave function at high momentum, we ignore the Pauli exclusion between 
different closed-channel bound-states, 0o, while keeping Pauli exclusion between channels through the factor 
u^. The higher order correction in this replacement is probably comes in as nonliner relationship between 
/i2k and (/)ok or extra terms (pi in expansion of /i2k- However, as we will show, /12k is always much smaller than 
1. These higher order correction is relatively minor and only modifies the conclusion quantitatively without 
affecting the result qualitatively. Furthermore, this replacement also solves the renormalization problem 
automatically because the two-body wave function decays fast enough at high momentum combined with 
the momentum-dependent interactions and therefore we have no more divergence in integration. 

Eq. 15.421 can also be understood in terms of BCS-type anstaz Hk (^k + "kOk^Lk + ^^^'L'^^-v) |0) (Eq. 
15. 2p which connects to coherent states of the two-body wave function (Eq. [O]), (j)^"^'^ + It is not hard 

to show ft,2k — (flkC-k) — ^ik^k for this ansatz. Using the relation (fy^'^'^ — Wk/wk(Eq. 15. 5p . we have 

/i2k = WkiWk = Mk( — ) cx u^'/'k"''' (5.44) 



42 



In a two-body problem of the Feshbach resonance, we expect 0'^'"^) oc 0o (see Chapter [3]). And we have 



h2k oc 0okWk 



(5.45) 



just hke Eq. 

To the first order of C, (/)ok ^ 'jrT^fj^, where k is the momentum scale related to the binding energy or 
the tuning, k'^ /2m ~ Ejj ^ rj. In the many-body relevant region where momentum is not significantly larger 
than the Fermi momentum, this quantity is actually very small and approximately a constant, -ij, because 
kp. The part where the many-body factor significantly alters the wave function is actually rather 
small comparing to the whole spread of the wave function over the momentum space, in the order of n 
{El, = ^^K^/2m). Within the range where differs from 1 significantly, the wave function (/lok is very small. 
(See Appendix IB3|) 



This fact is valid across the full range of the crossover. At some regions of crossover, A2k (or a) may be large 
and the total closed-channel weight is comparable or even larger than that of the open-channel. However, for 
each energy level in low momentum, the factor 0ok is so small that the product of /)-2k is still much smaller 
than 1. In summary, the full closed- channel (including normalization) is always small in the low momentum 
region even when the total closed- channel weight is large. This justifies our following perturbative treatment 
on the inter-channel Pauli exclusion between two channels. 

Not far away from the resonance point, absolute detuning 77 is always close to binding energy E^'^ of the 
closed-channel bound-state, (/jq. Using Eq. 15.411 we can write 




A2k = /i2k 



(i^lk + -E^3k) 



/l2k 



(£^k + a + ?7) 



Combining the above equation with Eq. 15.401 we have 




We replace ft.2k with au^(/)k as in Eq. 15.421 The above equation becomes 




k' 



k' 




k' 



k' 



k' 
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The last term can be rewritten as V'pk'/'ok = (2ep + £^b)0op using the two-body Schrodinger equation of 
the isolated closed-channel Eq. 15.431 We move the last two terms from r.h.s to l.h.s, and using = fk ^ A* 

a^k ^-^^6 + ?? - 2/i + -Ek - Ck + ^ Vkk'Wk'^k'^ = ^ ^^kk'/iik' 



Multiply both sides with (j)^ and integrate over the momentum, 



a 



k kk' 



^Eb + l]-2^l + J2 't'liEi, - a)0k + J2 <^k^'k'^4k'0k' 

k 

We rewrite this equation in the form 



kk' 



^^Fkk'/iik' (5.46) 



Ekk'</'kykk'feik' 

Ai(??) = - X cl>l{E^ - a)</'k - <^kt'k'Vick'0k' (5.48) 

k kk' 

Comparing this to the two-body problem Eq. 13. 6[ detuning is shifted by the many-body effects /i (mostly 
due to the Pauli exclusion within the open-channel) and Ai (mostly due to the Pauli exclusion between 
channels). Ai depends on detuning rj (through E\^, Uk), but the dependence is rather weak because the 
integration is over mostly the short-range quantities and insensitive to the detuning. We will discuss more 
detail about it later in Sec 15.4.2] Using Eq. I5.47[ we can express the a in /i2k (Eq. I5.42p . and then plug it 
into Eq. [07I 

A X^TT h ^2 V-r^. Ekk'p'ypp'V'/'k'ykk'^^^fcik 
^ip = ^pk'^ik + 2^ >^p'«0p' V = ^pk^ik + (_^^^_^_2^- Ai) 

Comparing this with the two-body problem, we can see that the detuning part (denominator of the second 
term) is shifted by 2/i-|- Ai and there is an extra u^, term introduced as the many-body effect. Nevertheless, 
none of these affect the high-momentum behavior. Therefore, the equation can be renormalized exactly as 
in the single-channel problem by introducing the long- wave- length s-wave scattering length as- We rewrite 
the above equation 

. x^frr ^ Ek'p'ypp'V'/'k'ykk' ^^ 

_ Ekk'p' ^PP' V'^k'^kk't^p.felk 

{-Eb + r|-2^l-\l) 

The second term in the r.h.s. has no divergence at high momentum in 3D due to the extra v^, factor. Actually 
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the factor, v^,, decreases quickly over a small range in the order of "gap" Aip; therefore, the summation 
in the second term is essentially only over low-momentum, and is very small. In addition, considering the 
short-range nature for Yppi, this term varies slowly over momentum, p. 

Multiply both side with (1 + TG), where T is the scattering matrix for the open-channel, and G = 
{(jj — Hq)^^ is the Green's function for a free pair in the open-channeiy 



(1 + rG)Ai = -Thi - Aa 

A.P . + TGU ^--'-:^^^^^^^^^^ = " Ed + TGUY,p,c^p,vl, (5.49) 

p pp' 

In principle, Aap depends on momentum, however, like Aip, it is approximately constant at low momen- 
tum. We are only interested in the low momentum/frequency properties as the system is cold, dilute and 
governed by only short-range interactions. The scattering matrix T is approximately its s-wave zero-energy 
value Anaslp,, Ai)/m although ds{fJ-, Ai) is the s-wave scattering length of the shifted (by 2/i -I- Ai) detuning. 
The s-wave scattering length, as, in the two-body problem without many-body shift 2/i + Ai is simply Eq. 

mm 



S ^ 

Here in the many-body context with the extra detuning 2/i -f- Ai, is 

a. - + ,_,^_,^ ) (5.50) 

where K, is introduced in Eq. 13.121 of Chapter [3] The Green's function is simply l/2ek at zero energy. 
Similarly, Aip and Aap is approximately constant at low momentum. We drop the subscript p of them, and 
pull Ai out of summation. 

Ai = -TYih, + GAi) - A, = _4vra.(/.,Ai) y ^1^ + ^ + ^ ^ - ^) - A, (5.51) 

Here the last step use the zero energy value of the zero-energy free pair Green's function G{lu = 0) = 
(— 2ek)~^. Dividing both sides with Ai, we have the renormalized equation 



1 = - 



m ^n^ik + ^2k)(£^ik + S3k) 2ek' 



Ai 



(5.52) 



^Here we use the relation of the scattering T-matrix, the free pair Green's function G = (cj — Hq) ^ and the bare interaction 
V = -U.S. 

T = V + TGV = -(1 + TG)U.fi 
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Note that as(/i, Ai) corresponds to the two-body s-wave scattering length at detuning shifted by 2/i + Ai. 

Now we can expand the first term in the parentheses to the first order of C, using Eq. 15.281 Keeping in 
mind i?k ^ ?7 at low momentum where summation is about, we havq^ 



1 = - 



47r5s(M,Ai) yr^fj_ I . 

m ^^2E^ 2ek 4.{^k + v)El' 



^ (5.53) 



The correction term (the third term in parentheses) does not have divergence in summation of high- 
momentum. In summary, there are several difference of gap equation here comparing to single-channel 
problems: 

1. The shift of 2/i in detuning through ds', 

2. The extra shift of Ai (see Eq. I5.48|) in detuning through ds', 

3. The extra term in Eq. [5331 

4. The extra term ~4{^^^f^fE^ ~ — ^^^j C in the summation of Eq. I5.53( 

Item [1] corresponds to a many-body effect common to both three and four species problem. In a many-body 
system, most of the properties can obtained when considering what happens at the next extra particle or 
particle pair. In a two-body system, we can pretend that the vacuum is the initial state and the next particle 
pair comes in from energy zero. In a many-body system, this starting point is no longer the absolutely zero, 
but the chemical potential. Excitations are counted from the chemical potential, fi. The energy detuning in 
Feshbach resonance is no exception. In a many-body formula, the detuning is shifted by 2/i (for a pair). In the 
broad resonance, the Fermi energy is smaller comparing to the detuning during the relevant probing region. 
Consequently, this shift is negligible. In the narrow resonance, the Fermi energy, however, is substantial 
and a detuning of the chemical potential needs to be taken into account. At the BCS side, /i is positive 
and close to the Fermi energy. This shifting just confirms the fact that most interesting phenomenon in 
a fermion system happens around Fermi surface instead of absolute zero. With an extra positive shift /i, 
detuning is reduced; therefore the crossover region is actually reached earlier (i.e. by larger detuning) than 
in two-body cases. On the other hand, moving toward to the resonance point and BEC side, the chemical 
potential becomes smaller and finally flips sign to negative. Accordingly, the detuning is reduced lesser and 
lesser; finally it is enhanced. Interestingly, the crossing point, fi — 0, is not shifted from a model with simple 
s-wave scattering length directly from two-body physics if only considering this shift (See Fig. 15. 2p . This 
effect is also studied extensively previously. (Sec. 6.2 of [Gurarie and Radzihovsky, 2007| ). 



■"Here we used ~ ~ W' 
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The next three corrections are unique for three-species problem where the inter-channel Pauli exclusion 
cannot be neglected. Furthermore, close look into Ai and A2 (Sec. I5.4.2[) reveals that they vary slowly 
across the full region of crossover and are functions of the density of atoms in the open-channel . They 
describe fundamental many-body effects, and do not have counterparts in a two-body problem. Both Ai 
(see Eq. I5.48P and A2 (see Eq. I5.49P involve overlap integrals between the open-channel wave function and 
the closed-channel wave function, (factor or describes mostly the open-channel wave function; while 
4> describes the closed-channel wave function). The larger the overlap of the two, the larger Ai and A2 
are. This has a very intuitive interpretation: more overlap leads more severe inter-channel Pauli exclusion, 
which in turn leads to larger Ai and A2 that describe this effect between two channels. In our model, an 
open-channel wave function is spread all over the real-space, (even at BEC-side, the real bound-state is very 
loosely bound comparing to the closed-channel bound state), while the closed-channel wave function is very 
sensitive to the binding energy, i?b(w ij). When the closed-channel bound-state is closer to the threshold, 
i.e., the binding energy is smaller, the closed-channel bound-state is more spread out in real space and has 
larger overlap with open-channel. Consequently, Ai and A2 are larger in such cases. Nevertheless, Ai is much 
smaller than the Fermi energy Ep, or the other shift, the chemical potential, 2/i. So the shift is not very 
large and it is still all right to treat it as a perturbation. The correction is shifted with the change of the 
density. In the many-body system, no dramatic jump in physical quantities happens at resonance due to 
the crossover nature, and it is hard to observe this shift over the resonance position though. 

5.4.2 Evaluation and estimation of Ai and A2 

From the last section, we see that a lot of the inter-channel Pauli exclusion is encapsulated in two parameters 
Ai and A2. Therefore, it is well- worthwhile to study them in more details. Ai is defined as 




k 



kk' 



Use relationship 



we can rewrite the above equation into 




(5.54) 



k 



k' 
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We can then replace the second term in the parentheses using the two-body Schrodinger equation of the 
isolated closed-channel Eq. 15.431 J2k ^ic^k = {Et + 2ep)0p. 



k k 

In the above summation, is only non-zero below or not much higher than the Fermi momentum. In 
this range, £^k,ek ^ V and we can neglect + Ck comparing to rj. Furthermore, 4>k varies slowly in this 
range, because the closed-channel bound state is much smaller than the interparticle distance. It is a good 
approximation to replace 0k with 0k=o ~ and to take it out of the summation. Put all these together, 
we can estimate Ai 

Ai « -T] |0fc=o|' ^ «k = -'7 \^k=a I' No (5.55) 

where No is the total number of atoms in the open-channel. It is easy to see that | Ai | is much smaller than i] as 
t=o|^ No < 1. Furthermore, using (j>k=o « \f^^ (Eq. [BJ3l of AppendixES]), we get Ai - '^Ep < Ep, 



\(Pk=o\ <^ i. rurinermore, usmg (pk=o ~ y 
where oq is the average interparticle distance and flc is the size of the close-channel bound-state. We have 
showed that Ai is indeed a minor correction over shift /i. It is not easy to estimate Ai precisely from 
microscopic parameters above because these parameters themselves are in turn hard to estimate precisely. 
Nevertheless, a key observation is that Ai depends on the (open-channel) density (ng = Vq"^ X^k^k) linearly 



to the lowest order of — , 



Ai « \^^\o (5.56) 



(0) 



We can perform the experiments at different densities and estimate the coefficient \\ 
A2 is defined as 

A2P . ;^(i + TG)pp ^Yl'f'-1''^"Zr'^ - " Ed + rOppFpp, v.^, s 

Here the argument goes similar as in the Ai case, considering the short-range nature, both Y^p' and 0p' 
vary slowly and can be approximated by their value at p' = for the momentum below or around the Fermi 
momentum. After pulling these two quantities out of the summation over p', we only need to sum over Wp, , 
which gives a familiar factor, the (open-channel) density, Uq. In addition, a — ^JTr^. The factor (1 + TG)pp 
integrate out the high momentum component and A2p is close to a constant at low momentum p. We have 

A2 ~ x'^^^Uoy^ (5.57) 
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In experiments, Ai and A2 can be measured at different total densities. The two coefficients, X{ and , 
wliich do not depend on the density, can then be estimated according to Eqs. 15.561 and 15.571 

5.4.3 Number equations 

There is one number equation for each channel. 



J2 G22e(-^""^-) = N, 



open 



J2 Gaae^-*'^"'-) = Ndose 

Note that the Matsubara summation formally diverges and we need to put in a small negative part into the 
summation in order to prevent the divergence. Here we put in a small negative part instead of a positive 
part because ^'2 and are conjugate of real particles, ^'2 = f/'b, ^3 = f/'c- The Matsubara summation can 
be performed with the normal trick of multiplying a Fermi function to summand and deform the contour 
(see Appendix IB. 21 for more details). For the summation at zero temperature, we just need to consider the 
positive roots, i^ik- It is straightforward to find 



? (^'Ik + i?2k)(£'ik + i?3k) ^ ^ 

AT ^ (-Elk - ^k)(-Elk +^k) - ^ ^Ik - ^1 (rt:q\ 

closed 2^ (£;ik + i?2k)(i?lk + i;3k) ^ (i?lk + i?2k)(£;ik + i?3k) ^ ' 



(£^ik-ek)(-Bik + ek + 77)- Ai 



Let us look at the equation of the closed-channel first, if we expand Ei k using Eqs. 15. 28115. 301 the lowest 
order is 

Here 71 is the correction for the fermionic correlation spectrum due to the inte-channel Pauli exclusion 
(Eq. I5.28p . We take the full value j^^^-^i instead of the low- momentum value because the summation has 
substantial contribution from high momentum. This is consistent with Eq. 15.381 and 15.411 if we assumqf 
A^cioscd ~ X]^2- From appendix IB.5| we know that the summand is much smaller than 1 in all regions. 
Particular, in the low momentum (< kp), the summand is (^^u^ ^ 1. Nevertheless, the weight spread in 
a very large range of momentum (~ rj) , so the resulting sum can be still in the order of total number N at 
certain detuning. 



*As discussed before, in the closed-channel /12k ^ 1 over all momentum. Furthermore, the summation goes over very large 
momentum range, the difference between /i2k and the closed-channel particle number is only large below or around Fermi 
energy, but the total weight of the closed-channel atoms in such a region is very small. 
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Note that the connection between A2 and the "Contact" in the theory of universahty ( [Tan, 2008a[ 
Tan, 2008b| ). The summand as density of momentum is vahd at momentum up to the scale l/flc, and is 



mostly determined by two-body physics except an overall factor at "high momentum" (kp <C fc <C l/j'c)- 
Furthermore, the contribution from the very high momentum (fc ^ l/''c) to integral is small. All the 
deviation there can be absorbed as a small correction in A2. To zeroth order of C, iJjk ~ and ~ 1 for 
the most of the integral domain. This summation can be written as 

A^closcd « ^ tL ^ \ (^-^l) 

^ (?k + ??)(2Ck + ?7) 

This equation has only one unknown parameter A2. Therefore, this equation can be used to estimate A2 
from experiments. The estimation of A2 from such procedure automatically includes the correction from 
high momentum. 

The open-channel number equation can be expanded perturbatively as well. 



k 



Ek-^k , 7lk (-Ek - ^k)(7lk + 72k) (i^k-^k)73k A: 



2Ei, 2E^ 4El 2Ei,{(k + El, + Tj) 2Ek{^w + Ei, + tj) 



^k - Ck ^k - Ck Ai ^ ^1 ^ 

2£;k 2£;k 2(£;k + ek + r/)^ 4El^ 



3 n (5-62) 



All terms behave well in 3D and do not need any renormalization. It is like the number equation in the single 
channel with a few correction in the order of C = A^/t^Ai. The second term, except the factor — v^, 

looks like the summand in the closed-channel number equation, Eq. I5.60[ which can sum up in the order 
of N. Nevertheless, the factor, decreases quickly in high momentum, becomes close to zero for k ^ kp, 
and limits the summation to only low momentum. This makes the sum a small correction as we discussed 
before. Similarly, ^ in the third term also limits the summation to low momentum and ensures that the 
sum is a small correction. At low momentum, this term is actually more important than the second term. 



5.5 Discussion of the mean- field solution 

As discussed before, the correction of the narrow Feshbach resonance can be taken into account in two steps. 
First, omitting the inter-channel Pauli exclusion, we only consider the chemical potential /i in the shift and 
the extra counting of the closed-channel. Then in the second step, we can correct the previous result with 
quantities originated from the inter-channel Pauli-exclusion unique in the three-species problem. 
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In the first step, the gap equation and the (open-channel) number equation are simplified to 

2E^ 



1 = - 



^open = E (5-64) 



k 



Here we only consider the shift of the chemical potential 2/i in 5,^ (Eq. 15.501) . 

In the second equation, we assume atg <C oq i'^o is the average interparticle distance) and we only study 
situations close to resonance, where /C ^ 5 or /i ((5 is the detuning from the resonant point in two-body 
physics). From the two-body analysis of the Feshbach resonance, we know the numerator in the Eq. I5.65l is 
closely related to the characteristic scale 6c defined in Eq. 13.191 



We can rewrite the as with respect to (5c 



In the narrow resonance, the weight of the closed-channel becomes substantial quickly. We can calculate 
the amplitude of the closed-channel according to Eq. 15.471 

{-Eh + 77 - 2^ - Ai) 

y, as a short-range interaction, only picks up the short-range part of h\(r) and 0(r) in the above integral. 
The short-range part of h\ is assumed to be proportional to the short-range part of the two-body wave 
function. Its normalization, however, comes from many-body physics. If we assume the zeroth order form, 
h\ = Ai/2£'k, Zhang has derived the form of /ii(r) in the real space, (Eq. (25) in [Zhang and Leggett, 2008] ) 

mAi 1 - r/as 

Comparing this one with the Bethe-Peierls boundary condition Eq. 12.41 we find that is the normalization 
factor. And Ai is determined by the number equation for open-channel density ng. Similar as the two-body 
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case (Eq. I3.2ip . we can write 



|2 



Extracting the normalization factor of /ii, the integral in the numerator then is just the same as that in the 
two-body physics and can be rewritten using the expression of /C according to Eq. 13.121 



0k^kk'/iik' 



kk' 



matg 



(5.69) 



Here we use rrir — ^m. We further assume that the energy difference between the resonant point (os — > ±00) 
and the level crossing point, where the uncoupled closed-channel bound state level is at the threshold 
of the open-channel, /C, is in fact much larger than fi or 6c, i.e. |— i?f, + r;| ^ 2/i, Ai. So we have 
\—Ei, + 77 — 2/i — All « /C (note that here 77 is equivalent to S in the two-body Eq. I3.23p . We can then 
write 

" ' (5.70) 



From Eq. 13.191 above, we can replace ICabg with ^J^E^ftFlm using Eq. 13.191 After dividing both sides with 
the total density ntot — ( — p — ) we find the weight for the closed-channel 

« - (5.71) 



ntota 128 ^j^ft)^^^ 

Pc is proportional to . This form is consistent with the formula obtained by Zhang previously [Zhang and Leggett, 2009] . 
In Fig. 15.21 we plot the open-channel weight 

/3o = 1 - /3c (5.72) 

We can see an example of the narrow resonance {Sc — O.OOlEp"'^^) in the Fig. 15.21 We rescale the 
detuning with Eq 



Eo = ^J28ji^lm ■ kp = 2^5cEf (5.73) 

Ai saturates at about 0.8£^f- At the BCS limit, rio ~ ntot, Po ~ 1- Ai is exponentially small, there is 
only very small amount of atoms in the closed-channel. Moving toward the resonance and then the BEC 
side, Ai increases and the closed-channel density increases as well. At the universality point of the two-body 



^We choose this scale in order to make sure that hakp would match askp in the single channel discussion in Chapter |4] if 
there was no —2/i shift as in Eq. 15.651 in order to comparing with the single channel case as in Fig. 14.11 When &/Eo = 1, 
—dskp = 1 if there was no extra — 2/i shift. 
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Figure 5.2: The narrow Feshbach resonance w/o the inter-channel Pauh exclusion vs. the single- channel 
model 

We plot the chemical potential fi, the open-channel gap Ai and the open-channel relative 
weight (3o for the narrow resonance; and chemical potential and the gap for the single-channel 
model. The gap in the open-channel Ai and the chemical potential /i are rescaled with the 
the Fermi energy of the total density, Ep°*\ The a;-axis is the detuning S rescaled with Eo 
(see Eq. I5.73P for the narrow resonance; while it is —l/askp for the single-channel curves. We 
have taken Sc = for the narrow resonance figure. We used the detuning from the 

resonant point for the a;-axis in the narrow resonance instead of —l/askp in the single-channel 
because the additional shift, 2^, considering in Eq. 15.651 Consequently, the chemical potential 
lines in both cases cross the a;-axis at the same point where ji = Q. 
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physics, Ai is in the order of the Fermi energy, and /3o is already substantially smaller than 1 if the resonance 
is narrow. Because the maximum of the ric is the total density ntot, we can see that Ai saturates. The 
narrower the resonance, the sooner it saturates. 

We can understand the saturation of Ai from another respect. Let us look at the open-channel gap 
equation Eq. 15.391 

Alp = C/pk/Mk + >pk/i2k 

k k 

Not far away from the resonance, the second term dominates the first term. And considering Eq. 15.421 
^2k = Q!0ok''^k' it '^ot hard to see that the closed-channel normalization factor a is proportional to Ai. 
Moving toward the BEC end from the BCS end, the closed-channel density increases as Ai increases and 
eventually dominates the open-channel density. In fact, when the closed-channel dominates, a w ^/N ; 
consequently Ai saturates and is directly proportional to the square root of the total density. 

Ai « VtV ^ rp=o, k0ok«k « V^t^iL (5.74) 

k 

Here the coefficient A^^^^^^ is common for one specific Feshbach resonance. It is possible to observe this 
saturation experimentally by measuring the fermionic excitation spectrum. 

We have the zeroth order result, /i'-^', A^"'' from the previous step. Now we can write down mean- 
field equations including the inter-channel Pauli exclusion effects, expand them and look for the first order 
corrections, ^^^^ and A^^-*. A2 is a first order quantity by itself, and is related to the closed-channel density 
by Eq. 15.601 Therefore, we do not need to include A2 when looking for the first order corrections. 

We start from the open-channel gap equation (Eq. 15.531) and number equation (Eq. I5.62p 

^ 1 1 X A^^k m A2 

47ra,(A*,Ai) ^^2E^ 2eJ ^ ^ + 7^)El 4™,(Ai, Ai) ^ ' 



2Ei,J ^ 2E^ 2(£;k + Ck + f7) ^ ^El 



We have rearranged these two equations to better compare them with the gap equation and the number 
equation used in the previous step, Eqs. 15.631 and 15.641 The main structure of the equations are the similar. 
In the gap equation (Eq. I5.53p . the second and the third term in the r.h.s are due to the inter-channel Pauli 
exclusion. In addition, in the l.h.s. also has a small correction due to the shift Ai(Eq. 15.661 and Eq. I5.50p . 

Ai (5.77) 



ai°)(^(o)) a.(M,Ai) yf25J9j^ 
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We can expand the l/(2£'k) with fi /i^*^) + fi^^^ and Ai — > A^"-* + A^^-*, where the zeroth order terms fi^^'^ 
and A^^'' satisfy gap equation in the previous step, Eq. 15.631 We then get the equation for /i^^^ and A^^'. 



= ek - (5.79) 



i^u-V^k + Ar^ (5.80) 

Similarly, we can expand the ft/ (2ii^k) in the number equation with ^ — > /i'"^ + /i^^^ and Ai — > A^^'' + A^^^ , 
where the zeroth order terms /i^'^^ and A^*^"* satisfy number equation in the previous step, Eq. 15.641 We then 
find the other equation for z^*^^' and A^^''. 

2^ -^k 2^ £;3 1 ^ 2i;k 2(Sk + a + »/) ^ AEl^ ^ ' 

Please refer to Eqs. 15.481 and 15.491 as well as Sec. 15.4.21 for Ai and A2. These two quantities describe the 
overlapping between two channels, i.e. the inter-channel Pauli exclusion. Their values are hard to estimate 
without the detail knowledge of the potential and wave-functions. Nevertheless, they can be derived from 
the experiments (Sec. I5.4.2p . Once we obtain them, we can solve the above equations for the first order 
correction /i*^^) and A^^' . In the above equations, all correction terms either have C explicitly or in the similar 
order. Given the non-singular nature of the crossover problem, we expect the correction /i*-^) and A^^' to be 
also in the order of C. 

There are three distinct length scales in the problem. The range of the potential, Tc, is the smallest. The 
intra-channel and inter-channel potentials are essentially zero outside this range. Potential energy dominate 
kinetic energy in this range and the correlation (wave function) is totally governed by the potential and 
might have large oscillation. The two-body correlation follows the two-body wave function in this range. 
The second range is the size of the closed-channel bound-state, Uc- The closed-channel bound state has 
only negligible weight outside this range. In this range, the many-body correlation still follow the two-body 
wave function in the closed-channel. We assume that the closed-channel bound state is close to threshold 
and therefore has most weight in this region. The inter-channel Pauli exclusion are mostly accounted in 
this region. Not surprisingly, this region contributes most in the integral of Ai and A2. When much larger 
than flc, many-body effects are important. However, the closed-channel has only negligible weight in it and 
therefore only the open-channel needs to be considered. 
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Chapter 6 

Excitation modes 



In the single channel crossover, fermionic modes are related to the two-body correlation of the original 
fermions; while the bosonic modes are related to the two-body correlation of the auxiliary new bosonic fields 
(order parameters (Ai, A2)). Similarly to the single-channel case, most modes are gapped with minimum 
at Ai and only one Goldstone mode for the in-phase phase fluctuation of order parameters (bosonic fields) 
is gapless with linear dispersion at low energy. 



6.1 Fermionic excitation modes and Bogoliubov transformation 



If we limit ourselves in the mean field level, we can interpret the transformation Tkik in Eq. (|5.2ip as the 
Bogoliubov canonical transformation, while the 3x3 fermionic correlation matrix B in Eqs. (I5.27ll5.30l) 
gives us the spectrum of the fermionic quasi-particle excitation. From Eq. (|5.14l) and Eq. (|5.2ip , the action 
is diagonal for new fermions (quasi-particles) field 



™.-k 
i'III,-k 



V'fi-k 

yV'c-k j 



Putting the above transformation into the operator language and using Aj, Xn, and Xm to denote the new 
quasi-particles, we have the relation 





-k 


<- 


-k 




-J 



= LiTl 



(6.1) 



56 



A1A2 



A1A2 



\ 







Uk 






Uk 










A2 











(E221) 



(lOTl) 



First of all, mixture of ak and foL^. (cL^.) indicates the pairing of atoms of the opposite momentum. Further- 
more, mixture of creators and annihilators dictates the approach of the grand-canonical ensemble, i.e., the 
ground state is a number non-conserved state. In addition, L matrix cannot be separated into two channels, 
which indicate the mixture of two channels. At the mean field level, both Ai and A2 are taken as real and 
constant. The new Hamiltonian at this level is 



H - /(Ai, A2) + £'ikA'j|_^^A'j + £'2 k-^"/! ^k-^ii -k + -^sk'^'in ^k-^iii-k 
And the spectrum is just like we calculated in Sec. 15. 3[ Eqs. 15. 28115.301 We quote here again 



(6.2) 



Elk ^Ek + 7ik ~ -Ek + uiC 
E2U = Eu + 72k « i^k - "kC 
^^3k = ?k + J? + 73k ~ Ck + 77 



(lOOl) 



Here 71, 72 and 73 are correction due to the inter-channel Pauli exclusion. 

This clearly shows that X^y^, j^., and X^^^ are fermionic quasi-particle excitation modes with spectrum 
Ei k and l'^T^ is the Bogoliubov canonical transformation to transfer the normal fermionic modes into these 
elementary quasi-particle modes . Here we see in the excitation, different species of opposite momentum, 
(a, 6) and (a, c), mixed together to form the elementary excitation due to the paring in the ground state. 
From Eqs. (|5. 28115. 3U)) . we see that the fermionic excitation modes basically follow the pattern in the broad 
resonance. In the broad resonance where the closed-channel only modifies the interaction in the open- 
channel, there are three fermionic quasi-particle modes: two (degenerate) Bogoliubov quasi-particle modes 



in the open-channel, £'k = y^^ + Af as in ECS theory (gapped at Ai in the BCS-likc states, /i > and 



•\//i^ + Af in the BEC-like states, fi < 0) ; and one high fermionic excitation mode in the closed-channel, 
Ck + normal gas. In the narrow resonance, first of all, the mean-field value of the gap Ai (or 
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\/ l-fi + A^) itself is modified as described in the last chapter (Sec. 15. 5p . Once that is taken into account, 
the above conclusion is approximately correct except high-order corrections in ^. The originally double 
degenerate excitation modes, i?k — vCk+^i ' '^'^^ sp^ii by C^i; while the third high excitation mode 
corresponding to the normal fermionic excitation in the closed-channel has a small correction in the same 
order. On the other respect, as discussed earlier in Sec. 15.31 corrections due to the inter-channel Pauli 
exclusion do not vanish when the inter-channel coupling, Y , approaches zero. 

6.2 Collective excitation modes 

Fermionic modes are derived from the correlation function Q^^ of the fermion fields '5, and therefore are 
mostly single (quasi)particle like. On the other hand, order parameters (Ai, A2) are defined in terms of 
collective behavior of many fermion atoms. Fluctuations of order parameters thus marked the collective 
excitation modes of the system. Here with a two-component order parameter, four independent modes exist: 
two for magnitude variation of each A^, internal phase between two A^, and the overall local phase 0(x) of 
Ai and A2. The first three change the magnitude of action and therefore massive; while the last one leaves 
the action invariant and thus massless. We study two modes of the phase fluctuation. The in-phase mode 
is the counterpart of the Anderson-Bogoliubov modes in the single channel problem; while the two channels 
introduce a new out-of-phase mode. 

For the single-channel crossover. Sec. 14.21 considers all (magnitude and phase) modes of the order pa- 
rameter fluctuation at the same time, and only calculates the low sound-like part of the spectrum. In the 
two-channel case, a general analysis with all modes in becomes unwieldy. We instead focus on one mode a 
time. We isolate the change in one mode and leave others at zero. 

6.2.1 The in-phase phase fluctuation 

The action of A, 5(Ai, AJ (Eqs. 15.191 I5.20p . is invariant if the phases of Aik and A2k rotate simulta- 
neously. We therefore conclude that there exists a massless (Goldstone) mode corresponding to the local 
phase invariance. We first study the massless two-channel in-phase phase fluctuation. Introduce the phase 
fluctuation 9, 

A,(x) ^ A,e'2e(.) ^ A,e-^2«(") 

This is equivalent to phase rotation of fermionic variable if) 
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Again, the phase 9{x) is common for both channels. With a phase shift, we can rewrite the action (taken 
mean- field value A(") = (Ai, Aa)^) (here we follow treatment of Nagaosa [Nagaosa, 1999| ) 



j 

+ A(o)t[/-iA(o) _ (^^)A(o) - A(°)t(^^)| 
i 



(6.3a) 

(6.3b) 
(6.3c) 

(6.3d) 



Note that here A^°) is the mean- field value of A^, a constant 2-component vector, no longer a functional 
variable. Here we see that 5*0 has the same form as before except it only takes the mean field value of A 
and it is described by the same correlation Gq (Eq. I5.19p . We can regard and 52 as perturbation for 
the so-called gradient expansion on V6' [Altland and Simons, 2010] . It is then obvious that 5*1 is in the first 
order while S2 is in the second order regarding with the (time or space) derivative of 9. Use the same spinor 
representation as before, S is bilinear of 4' a-nd therefore we can formally integrate out ip. 



S[e] = const. + In det^" 1(61) 



(6.4) 



We write out the formal Green's function according to the above action (with respect to the Nambu-like 
spinor) 



1 



k.k' 



( wi. y\ , 1 . (k-kO.(k + kO ^^,^ 



2m 



K2k^k' ^^y^^M- q')^(9)^('?')'5(9 + + - fc')^3 
2m pV 



(6.5a) 
(6.5b) 

(6.5c) 



where Gq is the same as (Eqs. 15.171 [5T9| . Here k = (a;„,k) and fc' = (aj„',k'). Like in the single channel, 1 
is 3 X 3 identity matrix, and 



/l 



CT3 




0-10 
0-1 



(6.6) 
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As the expansion in Eq. 14.201 we can look for the expansion Q ^ over Ki 2- 



trhig-i =trhiGo V tr(Go^) - ^tr(Go^Go^) + • • • (H:^ 
For the first order, tr(Go-^)j 

tTiGoKi)=Y,GokKikM=0 (6.7) 
fc 

tr(Goif2) =J2GokK2k,k 

k 

= - ^ ^ E ^'^(^o '^'^s) E 

2m pV ^-^ ^-^ 

k q 

= -£^E'^'%)^(-'^) (6-8) 
q 

Here we use the fact tr(Go fcCs) — n. tr(Go-fir2) is already in the second order of 6, and we only need 

to keep the expansion of K2 to this order. On the other hand, we have to go to the second order of Ki for 
the second order of 9. 



tr(GoXiGoXi) -Etr(Go .k^qKi k^q^kGo kKi k.k-\-q) 
k,q 

^ E 6{q)e{-q)[{-ujl)ir{Go.k+,a^Gok(Jz) 

fe,g=(cj,„,q) 

^ E + |-)(fcj + ^)tr(Go,fc+«G'ofc) 



i,3 = {x,y,z) 



(6.9) 



(6.10) 
(6.11) 



Here we introduce q — (a;„i,q) — k — k' . We can take q — in n^^'^lq) and 7r(^)(g) for low frequency and 
momentum (oJm, 2^IqP ^ ^1,2/ Use the previous result of Sec. 15.31 we can calculate the Green's function 
in the lowest and the first order of C, in tt's. After some long but straightforward algebra (see Appendix. 
lR4l). we find 



(6.12) 



7r(^)(0) = 



(6.13) 
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Combining all these together, we have a new action for the phase fluctuation 9 



(6.14) 



The correlation determines the velocity of the Anderson-Bogoliubov collective mode. The second term in 



Eq. (EUD, -E, 



2£S(g^+i)) ~ ~ Sk ^'^^y correction in the next order. Thus, this mode 

is qualitatively similar as the one in the single- channel with some correction in the order of C <JC 1 (see 
Appendix |B.5|) . 



6.2.2 The out-of-phase phase fluctuation 

Another interesting phase fluctuation of order parameters is the phase fluctuation in two channels out of 
sync. This mode is unique for a two-channel system without any counterpart in the single-channel model. 
When the phase fluctuation of two channels are out of sync, the inter-channel coupling strength changes. 
Thus, it is expected to be a gapped (massive) mode. We narrow down to the mode that the phases of two 
atoms (V'fc and i/'c) are opposite and leave all other modes constant. 



ipbix) 



?/;fc(a;)e+*'*("=) 



Ipbix) 

\Mx) J 



^ Mx) ^ 

^h(a;)e~^^(^) 
yV;c(a;)e+*^(^) J 



The order parameters do not have a simple transformation because they are connected to two channels via 
2x2 interaction matrix U, which mixes two channels (Eq. 15. lip . 




This term cannot be easily written in terms of mean- field value Aj. On the other hand, as mentioned before, 
we freeze all other modes to their mean- field value except 9. We therefore use another two-component of 
atom pairs(V'6?Aa, '4'c4'a), which is the linear recombination of (Ai, A2). 

Comparing to the in-phase phase rotation, the out-of-phase one does modify the action's magnitude. 
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Write the interaction term in the momentum space, 




-i0 



I)hlh2 + Y*{e 



l)hih2 



Note that here we take the mean-field expectation of pairs of fermions hi — (i/'bV'a) = X^k^ik; ^2 = 
(V'cV'a) — X)k ^2k, and their complex conjugatefQ. Comparing with the action of the in-phase phase fluctu- 
ation (Eq. 16. 3p . we find one more terms in the new action. 



j=(a,b,c) 

Si[9,-4>i,4ii] = / dx^i{drO){-4jh%IJh-4>c4>c) 



V0 



Zmi J 



s^m ^ dx Y{ 



-i2e 



I)hlh2 + Y*{e+'^^ ~l)hih\ 



(6.15a) 

(6.15b) 
(6.15c) 
(6.15d) 
(6.15e) 
(6.15f) 



Let us look at 6*3 first. We can expand it around small 9. 



53 [61] = / dx Vie-'''" - l)hlh2 + r*(e+"'' - l)hih* 



dx^i2 Yhlh2~Y*hih*2 6-2 Yhlh2+Y*hih*2 



From Eq. [QTl 



hik = A 



{Elk + -E2k)(-£'ik + £^3k) 



(6.16) 
(6.17) 



If Ai — >■ Aie"'', all /iik has argument (p (or Tr + (p). Thus hi has argument (p (or 7T + (p) as well. Also consider 



^■i/)-!/) and ipip are complex conjugate in the mean-field level. 
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the mean- field equation Ai = Uhi + Yh2, We can rewrite 

Yhlh ^ hlAi ~ hlUhi (6.18) 

U is real for a Hermite interaction. So this quantity is real. In addition, Y*hih2 is complex conjugate of 
Yh2h\. So they are equal because they are both real. From all these, we conclude that the linear term of 
6 vanishes. This is actually expected for a perturbation around the saddle point. Now we can rewrite this 
term 

S:i[6] = I dx- AYhlh2e'^ ^ 0{q)e{-q){-AYhlh2) (6.19) 
•' q 

Yh\h2 is actually the expectation of coupling for the two channel correlation. It is not hard to see that this 
value is negative in the minimum (saddle point). 

The other two terms in the expansion of the action, Si and ^2, go in the same way as the in-phase phase 
fluctuation except an extra i because only hyperfine species (h and c) participates. They actually produce 
the same result except only a half of the particles (6, c) participates, (n^ -\- ric ~ Ua ~ n/2) The full action 
about 9 is 

S[e]= J dxY,e{q)ei-q)[\n^''\0)iu;l-ul)^£-q'] (6.20) 



where 

2 _ mYhlh2 



(6.21) 



7r(o)(0) 

As expected, the out-of-phase phase mode is a gapped mode, with a spectrum starting from ujq. When two 
channels' phases fluctuate out-of-phase, the original self-consistent mean-field equation cannot be satisfied 
and this fiuctuation is coupled with the fermionic modes. Its general spectrum is expected to related to the 
pair-breaking energy, but Eq. 16.201 can only be calculated numerically in general. Particularly, Yhlh2 is 
hard to estimate except at the ends of crossover. At the BCS end, Yh\h2 is close the expectation of the 
interaction energy calculated according to the effective interaction in the open-channel. We can estimate it 
as Yhlh2 ^ - iV(0)Af, where iV(0) is the density of the state at the Fermi energy. 7r('')(0) A^(0). So we 
can estimate that ojq ~ Ai. This mode is inside the fermionic excitation spectrum. At the BEC end, the 
interaction energy in the open-channel can be estimated roughly as A'opon \lA\ 7r*-°-'(0) ^ VqA^/ Note 
that here the open-channel atom number iVopen might be significantly smaller than the total atom number. 
Here, Ai ^ nV^as and ~ . Put all these together, luq ^ n^^^ ImI^^^ /Ai ~ a^^ . It is in the order 
of the ionization (pair breaking) energy again and is only related to two-body physics. 
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Chapter 7 

Conclusions 



In this thesis, we study the narrow Feshbach resonance in the three-species case where two channels share 
the same species. 

In general, for the narrow resonance without a shared species, the main correction to the single-channel 
result comes from the extra counting of the atoms in the open-channel, which leads to the extra shift 2/i in 
Ug, and the closed-channel, which leads to the extra number equation. Two number equations exist, one for 
the open-channel and one for the closed-channel. The open-channel number equation resembles the number 
equation of the single-channel model. 

When there is a common species, however, the Pauli exclusion between two channels due to the common 
species in the three-species narrow resonance, calls for careful consideration. Our treatment follows the 
spirit of the so-called "universality" idea. For a dilute system with a short-range potential, such as the 
dilute ultracold alkali gas, the short-range part of a two-body correlation does not significantly change from 
two-body to many-body. This particular feature justifies using the two-body quantities (e.g. the s-wave 
scattering length a^) as the boundary condition for the many-body correlations. The most part of the 
correlations, where no two particles are close in short-range, is essentially free from interaction. This long- 
range part of the correlations does change from two-body to many-body. Unique in our problem is the 
involvement of the Feshbach resonance. A resonance makes the system extremely sensitive to even small 
change in the relative weight between two channels. We use the fact that, in the Feshbach resonance, the 
closed-channel bound state, as a short-range object, is still not very sensitive to the resonance; hence, 
we can, within the universality idea, apply a simple boundary condition to the closed-channel correlation. 
However, the indirect interaction mediated by the closed-channel dominates the direct interaction within 
the open-channel. Consequently, a small change in the closed-channel bound-state weight does affect the 
short-range wave function in the open-channel, originally expressed through the boundary condition using 
the s-wave scattering length Ug. 

When the spatial extension of the closed-channel bound state is of the order of the inter-particle distance, 
ao, or even larger, the Feshbach resonance in the many-body context is a genuine three-species many-body 
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problem and no simple solution is available. By contrast, when the bound-state's spatial extension is much 
smaller than the interparticle distance, the ratio (flc/ao) serves as the expansion parameter and we can extract 
the effect of the inter-channel Pauli exclusion perturbatively. In essence, we can then ignore the many-body 
effects within closed-channel bound-states, while only taking into consideration the Pauli exclusion between 
channels. A few new parameters need to be introduced and can be calibrated from experiments, such as 
Ai (Eq. 15.481) . A2 (Eq. I5.49p . Mean field properties can still be determined through gap equations and 
number equations similar as in the single-channel case. The excitation modes are also close to the original 
single-channel result with correction of the order of Oc/ao, where Oc is the spatial extension of the uncoupled 
closed-channel bound state and aq is the average interparticle distance. 

We can distinguish three different regimes for the many-body energy scale Ep when compared with the 
two-body physics. 

1. The regime where Ep is much smaller than the characteristic energy 6c (See Eq. I3.19l in Sec. [3]). Around 
resonance, the closed-channel weight is negligible in this regime. We can then ignore the closed-channel 
completely and approximate the effective open-channel interaction by a pseudo-potential characterized 
by the s-wave scattering length, Ug. 

2. The regime where Ep is larger than 6c, but smaller than the closed-channel bound state binding 
energy Eh, or the bare detuning r/. The closed-channel weight is significant in this regime and cannot 
be ignored, which means Jiopcn + "-close = "-total- Also, we need to take into account the fact that, in 
the resonance formula of a^, shifting should be counted from the Fermi surface instead of from zero 
as in two-body physics. These effects, considered previously [Gurarie and Radzihovsky, 2007| , are also 
carefully analyzed in this thesis. Furthermore, the Pauli exclusion between two channels needs to be 
taken into account when these two channels share a common species. Nevertheless, at low momentum, 
where the open-channel has most of its weight, the closed-channel bound state only has very small 
weight. Consequently, we can still treat the Pauli exclusion between channels perturbatively using the 
expansion coefficient Ep/rj. 

3. The regime where Ep is larger than the binding energy E^ or the bare detuning 77. We have a genuine 
three-species many-body problem. This regime can be achieved in two ways. One way requires a very 
large Ep, which indicates a very dense system. However, it would be hard to imagine that the original 
dilute alkali gas model still applies in this case. Various approximations in the model would probably 
break down beforehand, such as the pseudo-potential approximation with as. The second way requires 
a very small ry, which indicates a genuine three-species many-body problem. This remains an open 
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problem. 

This thesis focuses on the second regime. Careful analysis leads us to distinguish two related but dif- 
ferent concepts, the total weight of the closed-channel and the weight of the closed-channel in a particular 
momentum level. Within our assumptions, namely, the spatial extension of a closed-channel bound state is 
much smaller than the interparticle distance, a very low occupation in low-momentum levels (fc < kp) of the 
closed-channel persists despite the large total weight of the closed-channel. A suitable small parameter Ep/i] 
(or ao/flc), upon which a perturbation theory can be developed over a non-perturbative zeroth order solution, 
emerges with this discovery. Pairing in many-body problems is known to be non-perturbative and has to 
be handled with proper non-perturbative techniques. Nonetheless, we can concentrate the non-perturbative 
part into the zeroth order (broad resonance) solution, while treating the other corrections perturbatively. 
The resulting theory handles the two channels in steps self-consistentlj^. 

In summary, a two-channel model has a two-component order parameter (Ai, A2): one component for 
each channel. The order parameter for the closed-channel can be determined by the number equation of the 
closed-channel (Eq. 15.611) 

Where = h?k'^ /2m — /i and 77 is the energy difference between two channels. The renormalized gap 
equation is given in Eq. 15.531 

where E'k = \/ ^i- Here, as{^,, \i) is the two-body open-channel effective s-wave scattering length with 
an additional shifting 2/i + Ai as Eq. 15.501 

a. = «bg(l + J^^^) » 

Ai and A2 are two new parameters describing the overlapping of the two channels, and they can be calibrated 
from experiments (see Eqs. 15.48115.491 and Sec. I5.4.2p . The open-channel number equation is Eq. 15.621 



TV w V 

^ * open ~ / ^ 



3 1 



(1532)) 



Here C, — A2/Ai?7 <C i?k appears in multiple places as the small expansion parameter. The open-channel 



^Here by "self-consistent", we refer to the self-consistency in treating the close-channel according to zeroth order in the 
open-channel pairing. The BCS-type treatment of pairing in the open-channel is known to be not self-consistent. 
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number equation and the associated gap equation need to be solved self-consistently to get the mean field 
result (Ai and /i). One of the most noteworthy new phenomena is probably that the open-channel gap Ai 
saturates in the BEG side of a narrow resonance. 

There are three fermionic excitation modes. Their spectrums with the first order correction due to the 
inter-channel Pauli exclusion are given by Eqs. 15. 28115. 3(}1 

-Elk ~Ei, + ulAiC (lEsa 
-B2k « - «kAiC (EMI 

With a two-component order parameter, the bosonic collective fluctuation modes are rich. We explored two 
modes about phase fluctuations. The two-component in-phase fluctuation is massless and the low-energy 
one. It is similar to the Anderson-Bogoliubov modes in the single-channel problem with a small correction 
in the order of C. The new out-of-phase fluctuation is gapped and the minimum excitation energy is in the 
order of the pair-breaking energy (Ai in the BCS-like states, ^ > and fi'^ + Al in the BEC-like states, 
fi<0). 

In our approach, we take the broad resonance result (or the single-channel crossover) as our zeroth order 
solution, upon which the expansion is performed. It is however known that the simple BCS-type pairing 
treatment is not adequate to quantitatively describe the whole BEC-BCS crossover region. Therefore the 
zeroth order solution used in this thesis (simple BCS type ansatz or saddle point) can be improved through 
further theoretical development. Nevertheless, we expect the perturbative approach used here to build 
the narrow resonance from the single-channel crossover result to be still valid then. Once the zeroth order 
solution (for a broad resonance or a single channel BEC-BCS crossover model) is patched over with whatever 
advancement, the correction of the narrow resonance in such a parameter regime, can still be obtained with 
a procedure similar as the one discussed in this thesis. 

The theory we have developed here is for zero-temperature only. This limit simplifies the calculations 
considerably. However, an extension to finite temperature along the same spirit should be possible. The two- 
component order parameters should persist at low-temperature. At higher temperature, these components 
are likely to decay at different temperatures. The open-channel order parameter, associated with a much 
lower energy scale (of the order of the Fermi energy or even lower), should turn to zero first. Then the system 
becomes a normal gas with a Feshbach resonance. From the above discussion, many-body corrections due 
to the narrow resonance (both the intra- and inter-channel Pauli exclusion) seem to be agnostic to whether 



67 



the system is in superfluid state or not. Thus, we expect that these many-body corrections can be carried 
out in a similar fashion. There should still be corrections over the detuning due to the chemical potential 
and the inter-channel Pauli exclusion (Ai) although a system is more likely to be a Fermi liquid in the BCS 
side and a gas of normal fermion-dimer-molecules on the BEC side in such a temperature. 
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Appendix A 



The variational approach using the 
BCS-type ansatz 

We can also formulate the problem using the BCS-type ansatz. It is not difficult to calculate the expectation 
of the free energy as well as the wave function that optimizes the free energy within this ansatz. The 
optimization process gives us gap equations, which determine the wave function with the constraint from 
number equations. We will see that this method yields the mean-field solution consistent with the one of 
the path-integral approach. However, this method is mean-field by nature and difficult to be extended for 
studying collective excitations. 

We rewrite the Hamiltonian Eq. 15.61 in momentum space, and restore the momentum dependence of the 
interaction coefficients. 

H=Y. ^k«k ak + ^k^k &k + ^kC^ Ck 
k k k 

- Y C^kk'flk ^'ik^-k'flk' - ^kk'flk ctkC-k'flk' 
kk' kk' 

- Y ^kk'ftJfotkC-k'ak' - Y ^kVa^'cik'^-kflk 

kk' kk' 

Here we only keep the zero center-of-mass momentum paring terms as in the original BCS work [Bardeen et al., 1957| 
We take the free atom at zero magnetic field as the zero energy. 

Ck = fc^/(2m) + ry\ {i = a,b,c) 

77* is the Zeeman energy of the i"^ atom at magnetic field B. The hermiticity of the Hamiltonian imposes 

C/k'k = C/^k', Vk'k-VkV (A.2) 
We then introduce the BCS-type ansatz 

l*> = n ("k + Vi^albl^ + w^ald^) |0) (A.3) 

k 



(A.l) 
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|0) is the particle vacuum state. We require |uk|^ + I'^kl^ + |w^k|^ = 1 for normalization. This ansatz is 
constructed like the original BCS ansatz for superconductivity. An alternative ansatz would be Jlk(''^k + 
?^kak^Lk)("k + "f^kflkC^k) l'^)' which is actually the same as Eq. IA.31 because the cross term vanishes due 
to the Pauli exclusion on the common species. As the original BCS ansatz, this ansatz does not have the 
fixed particle number. We will just require the expected value of number operators match the total particle 
number. For all interaction terms, we get two types of contributions, for example, 

(t/kk'flk^-k^-k'ak') = ^kk \vkf + C^kk'Wk'<'Wk< 

k k#k' 

The first term is the Hatree term and the second term is the more interesting pairing term. 
This gives the full free energy as 

F={H- fiN) 

= E(^k + C^)l«k|' + E(^k + e£)kkr 

- X! ^kk |wk|^ - Y t^kk'Wk'"k'WkWk 



|2 _ 

kT^k' 



Vkk l^fkl^ - Y ^kk'^k'^ik'^k^^k 
k k^ik' 



(A.4) 



X^'^kkWkWk- Y ^kk'^k'^k'^^k^k 
k k^^^k' 



where we have set 



^kkU'k^^k- 2^ ^kk'W^k"k«k'"k' 

k kT^k' 



ca _ a _ a _ b _ b _ c _ b 

^k — ^k Mi Sk — ^k M J ?k — ^k M 



Two chemical potentials are added to make sure the Ua — ni, + ric = \n. In principle, ^Z'' does not need to 
be equal to /i'' as there is no exchange or conversion between (a) and (&,c). In fact, the structure of the 
ansatz guarantees that Ua = ni, + ric- Therefore, we set — /i** for simplicity and drop the superscript on 
chemical potential hereinafter. We will also drop the Hartree terms because they are only related to the 
density and thus can be absorbed into the chemical potentials. In the second summation, we will ignore the 
fact that the summation only goes through k 7^ k' because the corrections due to this restriction lead to 
only the higher order terms. 
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A.l Exact gap equations and number equations 

We introduce two new parameters ft-ik and /i2k, which corresponds to the mean field values (equal time 
average) of the anomalous Green's functions, (in the same way as hi^2 in Sec l5.4p . 

ul + vl + wl = l (A.5) 
UkWk = /ilk (A.6) 
UkWk = /i2k (A.7) 



We can solve Wk, Wk, w^k in terms of /iik and /i2k (all parameters arc taken as realjj. One complication in 
solving above equations is that Wk, fk and Wk are all monotonic functions of k, while ft-ik or ft,2k is not in 
the BCS end. So, we need to be careful when taking the square root. We introduce sgn^, for such purpose. 



1 



= o l + sgnfeVl-4/i?k-4/i| 



2k 

„2 2/l|k _ / 1 . /l _ /1?,2 _AU2 



„„2 ^"2k _ ''2k 



1 - sgn^ \/ 1 - 4/i2^ - 4hl^^ (A., 



l + sgnfeTT^IT^f— 4^ 2{hl^ + hl^) 



2k 



where sgn^;, = 1 for all fc in BEG cases, and for large k in BGS cases. In BGS cases, sgn^. = — 1 when k 
is small. In the single-channel crossover, sgn^. ~ sgn(ei; — /i), and the turning point is at = 1/2 which 
corresponds to zero chemical potential, /i = 0. The two-channel crossover is more delicate to treat. The 
turning point is still at = 1/2, but it is no longer exactly at /i = 0. The sgn^, function is very important 
in number equations Eqs. IA.141 They however can be absorbed later on, after we convert the formulas back 
into notations of (uk, Wk, w^) or S^k- 



^It is not hard to prove that the parameters that optimize the free energy are real, within an overall phase factor, when the 
interactions are real. (See also footnote |5] in Chapter [5] page 1391 ) 
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Eq. IA.8[ s derivatives over /iik arq; 



duj 2/iik 



2 

.2 



=sgnfe- 



The second equation is not very obvious, but can be obtained by noting that -gj^ — ~'gj{^ ~ Similarly, 
derivatives over /i2k are 

dul 2h 



k ^"•2k 

SgUfc 



^^'k ^ 8fe^i,fe2kSgnfc 



2 (A.IO) 



hlJi2V (l - sgUfcVl- 4/1^-4/1^1^)' 

2Vi-4/^^k-4;^^,(/.^, + /.y2 

_ 2/i2kSgnj, 8/i^i,/i2kSgnfc 



a/i2k Vl-4/^^k~4/i^, VI - 4/^^k - 4/^L (l + sgn. VI - 4/^L " 4/^^, 
We can obtain the gap equations by differentiating the free energy with respect to /iik and /i2k- 



/iikSgUfc b , 4/iik/iikSgnfc 

= Sk 



Vl^4)^?;:^4)^ ^1 _ 4/.^, - AKi^ (l + sgn, VI - 4/i^k " 4/i^k) 

- ^kk'^lk' - H i"kk'/l2k' = (A.lla) 



/i2kSgnfc 4/i2,^/i2kSgnfc 

= Sk 



Vl - 4/i^k - 4/i2k a/I - 4/i^k - 4/i2k (l + sgiifc Vl - 4/i' 



Ik 4/l2k 



- 5Z ^kk'/i2k' - 51 ^kk'/iik' = (A. lib) 

k' k' 

where 77 = — = rj'^ — rf is the bare Zeeman energy difference and is larger than most other energy scales, 



Here we follow the same convention by taking hik and /i2k as real 
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such as Ep. It is close to the binding energy of the closed-channel bound state when it is not too far away 
from the resonance point. We see one disadvantages of the ansatz method: unlike the path integral approach, 
where hi^ and /i2k have a very specific forms, Eqs. (|5.371 15.38p . and ready for further approximation and 
analysis, we do not have much guide for these two quantities here and need to make some guess. 
We can define order parameters 

Aik = J2 t^kk'/^lk' + J2 ^kk'/^2k' (A.12) 
k' k' 

A2k = ^kk'/i2k' + ^ ykk'/^lk' (A. 13) 

k' k' 

Both of them should have little k-dependence at low momentum for short-range interactions. We will drop 
the k subscripts in both A in the following. Two number equations can be obtained. One for each channel. 



^opon = 2^ «k = E fh2 ^ ( 1 - sg^fc VI - ^hf^ - ^hj^ ) (A.14a) 



2 



iV,jo.e = 2E^'k = E (h^ +\2 ^ ( 1 - sgHfeV^l - 4/1?^ - 4hlA (A.14b) 

^ V xk 2k/ \ / 

A. 2 Approximate solution of the mean field equations 

Here instead of solving these equations, we just demonstrate that they are consistent to the saddle point 
solution obtained in the path-integral method up to the first order of correction C. 

Let us look at the closed-channel gap equation (Eq. IA.llb|) first. It can be rewritten as 



/i2kSgnfc 



2k 



V 



(l + sgn,Vl-4/iL-4/iL) ) 



E^kk'/^2k' - E^kk'^lk' = 



We can put some terms back into the notation of (wk, Wk, ''^k)- 



sgnfe Jl - 4/1?^ - 4hl^ = {ul -vl~ wl) (A.15) 



l + sgn,Vl-4/i?k-4/iL = 2< (A.16) 

(A.17) 



iMk _!i 

,2 



l + sgn,Vl-4/i^k-4/iL 
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We rewrite the closed-channel gap equation using these relations 



''^ ' SCk + ^ - ) - ^ Fkk'/i2k' - ^kk'/iik' = 



k' 



At high momentum (3> kp), the third term in the parenthesis is negligible and the the denominator — 
— ~ 1. This equation is then very similar to the two-body Schrodinger equation (Eq. 13. 4p . For the 
closed-channel, the interaction, or the closed-channel bound state, 0O: has non-zero value for a range in the 
momentum space much larger than kp. Thus we expect the closed-channel correlation ft,2k oc (/)o/(k) where 
/(k) approach one at high momentum. This is in fact the same argument as Eq. 15.11 in the beginning of 
Chapter [5] The next thing to notice is that the third term in the parenthesis might dominate the other two 
terms in certain situation (low momentum, BCS) because the denominator can be very small in such 
regions. A nature way to compensate this big term is to sejfl /(k) = u^. Now we arrive the same conclusion 
for the closed-channel correlation as in the path-integral approach (Eq. I5.42p 

/i2k = a(/)okUk (A. 18) 

Note that a is not fixed yet and should not be identified as a in the path-integral approach (Eq. 15.421) a 
prior. The equation becomes, 

«^ok«^ - (2^^ + V-§v)-T, ^kk'5,/>ok4 - E ^kk'/iik' = 



k \ "k 



k' 



Using the similar approach as in the path-integral method, we multiply the equation with (j)Q^^ and integrate 
the equation in order to find a. In the denominator of the first term ^ t;^, wl. in the most of the 

integral domains (up to order of K{ri) in the momentum space). The equation is approximately 

- ^ok^ (2a + V) - ^kk'^/'ok^Sk^ - E '^Sk^lck'/^ik' = 



k ^ kk' kk' 



And it is not hard to find 



And Ai is simpler than Eq. 15.481 



Ekk' KYkk'hik' 

—Eh + 77 — 2^ — Ai 



'^liv) = - E <?^k^^k'^kk'</'k' 
kk' 



(A.19) 



in Chapter [5] is defined to be {^^ + i?k)/2i?k a^nd does not carry an obvious physical meaning as in this appendix. 
They are the same in the lowest order of our expansion over C, for the narrow resonance. This does not affect our analysis 
because the closed-channel correlation is actually a first order quantity. 
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But it should be equal to Ai from the path-integral approach at the first order of ( (Eq. I5.48P because this 
term is the dominated one (see Sec. 15.4. 2[ Eq. I5.55p . We will just use Ai hereinafter, a is the same as a in 
the path-integral approach if /iik here is the same as the one in the path-integral approach, which we will 
show in the following. Particularly for low momentum, (/)ok ~ TffeZ ^^^^ Appendix IB.3[) and it is not hard 



to find at low momentum 



a0ok = — A.20 



Now let us check the open-channel gap equation. We can rewrite it as 

/ 



/'-ikSgnfc 



2k 



^k' + ^^''^ -r, I = Ai 

y (l + sgn, - - 



Again, using Eq. IA.16[ we replace the denominator with 2m^. Also we replace /i2k with a^ok""^ using Eq. 



^ikSgUfc 



(2a + (5</'ok)'r/) = A 



At low momentum, the second term in the parenthesis can be simplified as A2/77 = C^i using Eq. IA.20I 



\/i - 4/iL - 4/i: 

We can solve /iik in terms of Ai. 



2 

2k 



, _ Af(l-4/^L) AKl-4/.^J ^ Al CkAi A, , 



(2^4- CAi)2+4A2 4£;2+4ekAiC 4^;^^ 



Here we take the low momentum value of /i2k ~ "''^kC the last step. This term is much smaller than 
because of the factor — at the low momentum. We can neglect this term. Now we have 



,2 Af ekAi 

'1"^ 4£;2(1 ^2 



(A.21) 



This is equivalent in the first order of C to the /iik in the path- integral approach (Eq. I5.37p . 

It is easy to see that h^^^ is the higher order as C and wy^ 1 all the time. So we always have ft-2k ~ Wk 
except in low momentum. However, the summation in the low momentum is very small, one order higher in 
C. The closed-channel number equation is then simply A'cioscd = ^k ~ Sk ^2k' '^hich is the same result 
in the path-integral approach as we discussed in Eq. 15.611 

We move to the open-channel number equation Eq. IA.14al At the momentum smaller than the charac- 
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teristic momentum of the closed-channel bound state k (ri ^ Eh factor "i^i . is approximately 

1, 

while on the high momentum 



A? 

A? 

< 1 



Therefore the factor "i ^ . limit the summation into low momentum. On the other hand, in the low 
momentum we have, ~ — thus we can neglect /i2k square root comparing to h^^^. And using 

Ik ^ 

Eq. ron 

So we have the open-channel number equation 

iVope„«5:i-|^-^C (A.22) 



This is consistent to the open-channel number equation derived from the path-integral approach (Eq. I5.62[ 
where the second term is negligible comparing to the third term). 

Thus we have shown that the mean- field solution of the path- integral approach can solve the gap equations 
and number equations of the ansatz method up to the first order of Q. 
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Appendix B 

Materials for Chapter 



B.l Diagonalization of the matrix Eq. (15.251) 



We need to find a unitary transformation L to diagonalize tlic matrix 



-Ek WfcAa 
+Ek VkA2 
\ukA2 UfcA2 +r]J 



(15:251) 



We drops ali tire k subscripts in tliis section because matrices in tlris section are decoupled in momentum 
and we only deal with one particular momentum k a time. We notice that the first term is proportional 
to an identity matrix and does not change by unitary transformation, we only need to concentrate on the 
second term. We rescale all elements with for simplicity in the following of this section. 



El, 



(-1 



R 






uy 


1 


vy 


vy 


t 



The secular equation of i? is {\xl — R\ =0) 

{x^ - l){x -t)- y^x + (u^ ~ v^)y^ = 



(B.l) 



We use 



1 here. We assume at the zeroth order, the three eigenvalues are —1, 1 and t. [t has 



weak dependency on energy as (^^ + ifj/Ek, however, at the low energy region of interest, we ignore ^fe.) 
Both y and t are larger than 1, however, we will verify that given condition y^ <C t, the correction is indeed 
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small and the expansion is reasonable (See Sec lB.Sp . Indeed, close-channel component can still be smaller 
than the open-channel component at low-k (in the order ofkp) due to the close-channel bound state is much 
smaller than the interparticle distance even when the total close-channel atom number is more than that of 
open-channel. And here all the quantities are about low-k unless specifically noticed. We expand the system 
to the first order of the dimensionless parameter = y'^/t (Eq. I5.26p . and find 



-1 
1 

t 



2'' 



Eigenvector 



uvy 
2t 



uvy 
2t 



uy 
t 



Try 
t 



uy yy 1 
t t 



Now it is easy to write down the corresponding diagonal matrix and the unitary transformation 

f-l ~u^C ^ 

l-v^C 



B = ioJnl + E 



\ 



(B.2) 



L 



uvy 
~2t~ 



\-- 



uvy 


uy 


2t 


t 


1 


vy 


t 


vy 


1 


t 



4cy 



uy ^ 



(B.3) 



Here L is not exactly unitary transformation, it is only unitary in the first order of C. We have 

B = iujnl + E L^RL + o(C) 

Alternatively, we can write L as 



L = I 





A1A2 

— U 



A1A2 



V 



A2 



(B.4) 



Here we use uv = Ai/2i?. 

In the above treatment, the small parameter C, is momentum dependent. If we restore the subscript k 
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and scale it back with i?k 



c 



^2 



(B.5) 



A momentum-dependent small parameter is not very convenient to work with, so we take its maximum value 
in low momentum (< Ep)- In the BCS-like states > 0), minE'fe — Ai, min^k ~ 0; in the BEC-like states 



(/i < 0), minEk = ■\/A^ + /? and min^k — We take the smaller values and have our expanding small 
parameter ((Eq. I5.26P 



C, = max C = 



^2 

Air, 



(B.6) 



B.2 Derivation of the mean-field equations ( 15.361 ) 

We have fermion correlation as a 3 x 3 matrix (Eq. I5.19p . 



IbJn -ik Ai 

Ai iw„ + 

A2 iiOn 



A2 





' All 


A12 


Ai3^ 


A* 
^12 


A22 










-433 J 



1 

JA\ 



Here we work in the momentum space, in which the system is nicely decoupled at least to the mean-field 
order. And we therefore drop all the k subscript in the rest of section. A general 3x3 matrix inverted as 
such, 

/ \ -1 / \ 

^22^33 —^12^33 —A13A22 

-AhA33 AiiAss-AisAh A^Ais (B.7) 

^-^^3^22 ^12^13 A11A22 - Ai2Ahj 

where \A\ is the determent of A. At the mean field level, all A^'s are real constants. We denote the mean 
field value of G as Go . The determent of Gq ^ can be expressed as 



|Go ^1 = {iLUn - Ei){iuJn + £^2)(jw„ + E3) 



(B. 
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where iJ^'s are defined in Eqs. 
can find the last term in 15.341 



And Go can be obtained according to the above rule. Now we 



tr 



/o 1 

0-00 

V 



^A*(iw„ + e + 77) 



{iUJn - Ei){ilxJn + E2){iUJn + £"3) 



(B.9) 



£1 + e + ?/ 



V A*- 



Here we perform the zero-temperature Matsubara summation in the third equal sign with the normal trick 
(see sec. 4.2.1 in [Altland and Simons, 2010 , sec. 25 in [Fetter and Walecka, 1971] , also refer to Footnote[2] 
at Page. [23l) . Because of zero-temperature, within three roots, —E2 and — i?3, we only need to take into 
account two negative roots —E2 and — i?3, assuming the correction is small. Similarly 



tr 



E2){E^+E,)-^ 



(B.IO) 



And we have 



((7-i)iiAi + (J7-i)2iA2 - ^ /ilk = 

k 

(?7-l)l2Ai + {U-%2A2 - ^ /i2k = 



Invert the interaction matrix U and we have Eq. 15.361 



B.3 The wave function for a short-range potential 

Here we discuss some possible generalization on the wave function for a short-range potential. This topic 



has been studied by Zhang Zhang and Leggett, 2009 . We will use some similar ideas. Outside the range Vc 



of a short-range potential, an atom is free and the Schrodinger equation is very simple. 

V^V = EiP (B.ll) 

2™ 

The equation has a simple solution for s-wave, — A'e^'^^ /r {k is imaginary for a scattering state). For a 
bound state, normalization A' is determined by connecting it with the short-range part of the wave function, 
ipo, and then requiring the full wave function normalized to 1. 
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Let us discuss the bound-state first, where k > 0. In the momentum space, there is also a universal 
behavior at low momentum, where krc <C 1. 



-nr 

V'fc = TT^v? / ^"-(^o + A'- )e-''^-'- 



,e 



(27r)3/2 

The first part for ip^ corresponds to the short-range part of the wave- function. 

The first term has very little k dependence for low memontum k <^ i/rc and the second terms is more 
important in this range. Furthermore, if the bound-state is close to threshold, the most weight is outside 
Tc, we can neglect the first term and we have universal behavior at low-momentum while the normalization 
factor A can be easily determined. 

Where Vo is the total volume of the system. And the wave-function is 



/ Sttk 1 / Sttk 1 , „^ 



The second approximation is when the momentum is low (< kp)- 

Besides the bound-state, if the interaction is weak and short-range, the low energy scattering state is 
well described by the s-wave scattering state ip cc 1/r — 1/a (Eq. 12. 4p , and its Fourier transformation in the 
momentum space has the similar form When considering many-body physics, in the low momentum 

below or around the Fermi momentum, wave function is modified by the many-body effect; but in the medium 
momentum, (still much smaller than l/^c), this l/Zc^ universal behavior is preserved. The distribution of 
particle in such momentum, fci? <C fc <C I/t'c, is 1/k*. This is actually the "high-momentum" (medium here) 
behavior {C/k^) described in Tan's work about universality [Tan, 2008a[|Tan, 2008b| . 

On the other hand, at very higher momentum (fc 3> l/?'c), the second term in the above is very small. This 
is because the smooth tail part of the wave function {(po) cancels out and contributes little in high frequency 
momentum oscillation. The high-frequency Fourier component in momentum space is solely determined by 
the wave function within the potential range (rc). This can be extend beyond the two-body wave function 
to the two-body correlation as long as the long-wave-length part is smooth. In all cases, two-body, or 
many-body, very high-frequency of two-body correlation follows the two-body wave function. 
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Incidentally, the two-body Schrodinger equation in the momentum space reads as 



2m 



(B.14) 



At very high momentum (determined by the interaction strength and potential range), the first term domi- 
nates, and we have the asymptotic form of the wave-function similar as Eq. IB. 131 



lim — A- 



1 



(B.15) 



1-2 2 

where — %^ = E. Note that this behavior is for a different reason and A is not necessarily equal A at low 
momentum discussed before. 



B.4 Evaluation of 7r(^)(0) and 7r^(0) 

Here we calculate 7r(°) and (Eqs. l6J0l KWj to the first order of C (Eq. using the expansion of the 

Green's function (Eqs. I5.33|) described in Sec. 15.31 



^(o)(0) =^tr(Gofc^73Gofe(T3) 

k 

« tT{T^B,'TlasT^B,'Tlas) 

k 

+ tr(TkSkB-^Tla3Ty,B-^Tla3 - T^B-H^Tla^^T^B-^Tla^ 



(B.16) 



+ T^B-'Tla,T^5^B-'Tla, - T^B-'Tla,T^B-H^Tla,) 
Note that k stands for both the momentum and the Matsubara frequency, (a;„, k). (5k is defined in Eq. 15.311 



5r = 



A1A2 



A1A2 


Uk 







Wk 


-Wk 






A2 

'7 



Introduce two matricies AI^ and M^., 



Mfe = Tla^T^B-^Tla^T^B-^ 
Mk ^ B^^Tla,T^B^^Tla,T^ 



(B.17) 
(B.18) 
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And we can rewrite tt^"^ (0) as 



(B.19) 



Here we use the cyclical property of the trace tT{AB) = tr(_Bj4). It is straightforward to calculate 

n \ 



/ a 



Ai 







Now it is easy to calculate the first term 

2A2 



+ 



(B.20) 



Only root ~E2k in the first term contributes in the Matsubara frequency summation at zero temperature. 



EMM.) = E 



2A2 



i?^ k + i?2 k) V V 2^^k (Ck + ?7) 



(B.2I) 



For the lowest order of the second term in Eq. (|B.16[) . we only need to take the lowest order of Bk 



iuJk-E]^ ^ 

iujk + Ei^ 
iouk + + V J 



(B.22) 



It is easy to verify at this approximation 



tr^^kMfe - (^kMfe) = 



(B.23) 



Combine Eq. (|B.21|) and Eq. (|B.23I) . we have 



. a2 



E' 



(B.24) 
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7r^(0) can actually be calculated exactly 

4(0) = Y.^h){k,MGokGok) (B.25) 

k 

tr(Go kGok) = J2 HT^L^B^^lItIt^L^B^^lItI) 

k 

= Y.HBk'B,') 

^ (B.26) 

=0 

B.5 Consistency of the expansion over C, 

In our treatment here, one crucial assumption in expansion is the smallness of A2/?/ comparing to 1. Here 
we check it. We have the closed-channel gap equation (Eq. I5.40p 

The first term on the right is relatively small comparing to the second term. We just keep the second term 
for estimation. Furthermore, we assume ft.2k = V^cipou., where (j)Qk is the normalized wave function of the 
isolated closed-channel potential satisfying Schrodinger equation (Eq. I5.43[) 

-i;^°Vop = 2ep0op-^F0ok (Eia) 

k 

Rearranging it, we have (especially at low momentum) 

^ V(po k = (2ep + Eb)(pop ^ ri<j>op 

k 

Here Ei, is the binding energy of the closed-channel bound state and 77 is the Zeeman energy difference. The 
second approximation is correct at low momentum (smaller or in the same order of the Fermi momentum) 
a.s Cp Eb Ri rj not too far away from the resonance. Put all these together, we have 

A2 « aEb(t)k=o 
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If we assume a simple exponentially decayed wave function as in Eq. EH from Sec. U 



1 /Sttk 1 ^ . 



Here Vo is the total volume and n is the characteristic momentum of the closed-channel bound state, 
rj ^ Eb = H^k'^ /2m. The second approximation above is only for low momentum. Collect all these together, 
we have 



3 ,3 

Sttk 1 jric I kpc\^ I Epc 



kpc is the Fermi momentum corresponding to density of atoms in the close-channel, which is much smaller 
than the characteristic momentum for the bound-state, k. Therefore we have A2 <C ??, even when Uc is close 
to the total density n. 

Now we check whether the corrections in Fermionic spectrum (Eqs. 15. 28115. 3D)) , are indeed small comparing 
to the zeroth order terms. 

Esk-ek + rj+^Ai 
Here, we mostly only concern of case of low momentum (k ^ kp). In Eq. I5.30[ 



2 



< 1 



^2 

Eq. 15.281 and Eq. 15.291 are slightly more complicated. Both of them involve ■p^^ at the BCS limit, the 
closed-channel density is small, fci?c is small and that makes this ratio small; when close to the (narrow) 
resonance, where ric is comparable to to the total density, at low energy, Ai is in the order of the Fermi 
energy, so does E'k- We have (we no longer distinguish kp c with kp) 

C=^-V^-^«l (B.28) 



More deeply, in the secular equation that leads to spectrum, Eq. IB. II We rewrite it without scaling to 
^'k, (We drop subscript k in the following equations for simplicity) 

{x^ - E^){x -^-ri)-Alx + AlE{v? - w^) = 
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It is not hard to use definition of u and v to find — v"^ = ^/E, and express E"^ = + A^. Therefore we 
have 

Here the first term is for free particles, and let us estimate the relative size of the last two terms. For 
low-momentum solution, we simply use Ai Ep, we find 

Al{x-i) ^AlEp^kp''^ 

This justifies our choice to neglect the last term when finding the lowest-order solution and then use the last 
term for correction. 

In another word, the above estimation is just saying that the total occupation number of the closed- 
channel at a low momentum level is much smaller than 1 in all regions of resonance (narrow or broad) 
because the closed-channel bound state is much smaller than the intcrparticle distance. This factor gives us 
a small factor, C ~ ~ ~ \/ "TTi upon which we can do the expansion. 
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